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I ■ Abstract 

O ■ Sr2Cu304Cl2 is an antiferromagnet consisting of weakly coupled CuO planes 

which comprise two weakly interacting antiferromagnetic subsystems, I and II, 
which order at respective temperatures T/ ~ 390K and T// ~ 40K. Except 
asymptotically near the ordering temperature, these systems are good represen- 
tations of the two-dimensional quantum spin 1/2 Heisenberg model. For T < Tjj 
there are four low-energy modes at zero wave vector, three of whose energies are 
dominated by quantum fluctuations. For T// < T < Tj there are two low energy 
modes. The mode with lower energy is dominated by quantum fluctuations. 



Our calculations of the energies of these modes (including dispersion for wave 
vectors perpendicular to the CuO planes) agree extremely well with the experi- 
mental results of inelastic neutron scattering (in the accompanying paper) and 
for modes in the sub meV range observed by electron spin resonance. The pa- 
rameters needed to describe quantum fluctuations are either calculated here or 
are taken from the literature. These results show that we have a reasonable 
qualitative understanding of the band structure of the lamellar cuprates needed 
to calculate the anisotropic exchange constants used here. 
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I. INTRODUCTION 

There has been a resurgence of interest in low-dimensional magnetism due in part to the 
desire to understand high-T^ superconductivity. The lamellar copper oxide systems, when 
suitably doped give rise to a family of superconductors with T^'s in the range about 30K.EI 
In these systems the Cu ions are essentially in a 3d^ configuration. Due to a large on-site 
Coulomb interaction, the states of this system which are accessible at ambient temperature 
have one hole per Cu ion, and hence the manifold of such accessible states is described by a 
spin 1/2 Hamiltonian having antiferrmagnetic interactions, which are strongest between nearest 
neighboring Cu ions in the Cu02 plane. That this system is a nearly perfect realization of the 
two dimensional (2D) spin 1/2 quantum Heisenberg model has been established by a wide 
variety of experiments.^ 

Recently, a variant of this system, Sr2Cu304Cl2 (2342), has been shown to display very 
interesting magnetic properties. ETO The structure of this systemEl is one in which an additional 
Cu ion (which we refer to as a Cull ion) is inserted at the center of alternate Cu plaquettes of 
the usual copper lattice, whose ions we refer to as Cul's. Although all the Cu ions are chemically 
equivalent, they play very different roles insofar as magnetism is concerned. The Cul's order at 
a relatively high temperature {Tj = 386K) and have properties similar to those of other lamellar 
cuprate antiferromagnets.0 With respect to the isotropic exchange interactions, the coupling 
between Cul and Cull ions is frustrated. As a result, the Cull's order independently at a much 
lower temperature, Tjj = 39. 6K into the magnetic structure shown in Fig. 0. For Tjj < T < Tj 
a very small residual anisotropic exchange interaction causes the Cull spins to have a small 
ferromagnetic moment, the study of whichQ led to the determination of the magnetic structure 
which has recently been confirmed by neutron diffraction.Q The study of the statics also led to 
the determination of several coupling constants in the Hamiltonian used to model this system. 

A natural continuation of this study was to investigate the dynamics of this system, and in 
the accompanying paperQ (which we refer to as I) an inelastic neutron scattering study of this 
system is reported. One interesting result of these experiments was that although the coupling 
between the Cul's and Cull's is frustrated in the mean-field sense, the spin-wave spectrum 
showed an incontrovertible signature of interactions between these subsystems. H'lI The nature 



of this coupling was described by Shender in a seminal paper.! Although this phenomenon has 
been identified in other materials,B the effect of this coupling, caused by quantum fluctuations, 
is perhaps the most dramatic in the system 2342, as described briefly previouslyEl and in more 
detail in I. As the Cull system orders for T < Tjj, the small gap spin-wave energies are found to 
increase sharply. This increase indicates that even though the CuI-CuII coupling is frustrated in 
the mean-fleld sense, quantum fluctuations lead to a signiflcant interaction between sublattices. 
A less obvious type of frustration arises with respect to the in-plane anisotropy associated with 
the bond anisotropy of the exchange interactions. When the moments lies in the easy plane, 
the exchange tensor for spins i and j in the plane has different values for directions parallel and 
perpendicular to the i-j bond. However, within mean fleld theory this anisotropy disappears 
when the average over all bonds is taken. But as before, there is a signiflcant residual interaction 
due to quantum fluctuations which gives rise to in-plane anisotropy. Finally, even classically 
frustration can be removed by exchange anisotropy which has a form similar to the dipolar 
interaction. We will refer to such exchange anisotropy as pseudo-dipolar. 

The purpose of the present paper is to calculate the spin-wave spectrum in order to give 
a theoretical interpretation to the data presented in I. From the discussion so far it is clear 
that most of these phenomena are outside the scope of linearized spin-wave theory. What is 
required is a nonlinear spin- wave analysis, i. e. an analysis which includes the effects of quantum 
fluctuations. In fact, from an analysis of the magnetic structure of the cupratescJ it was shown 
that there are several perturbations away from the linear analysis of the isotropic Heisenberg 
model that one must consider. These are the ones mentioned above, namely, a) quantum 
fluctuations of otherwise frustrated interactions, b) quantum fluctuations of the anisotropic 
in-plane exchange interactions, and c) pseudo-dipolar exchange anisotropy between the Cul 
and Cull subsystems. In a simplifled way, one can categorize these effects in the way they 
contribute to the spin-wave energies, which is given by the famous formulae^ 



u; = PHeHa , (1) 

where He {Ha) is the exchange (anisotropy) field and we work in units such that cj. He-, and 
Ha are all energies, usually given in meV (ImeV/fc^ = 11. 6K, lmeV/h= 241.8Ghz.) We will 
see that the out-of-plane anisotropy of the exchange interactions gives rise to a corresponding 



out-of-plane anisotropy field H'^^ which has been understood in terms of the out-of-plane 
anisotropy in the exchange interactions without reference to fiuctuations.ll3ll^ In contrast, the 
in-plane anisotropy of the exchange interactions, when summed over bonds, averages to zero 
and therefore only contributes when fluctuations are taken into account .liHlO The mechanism 
studied by ShenderS contributes to Ha except for the Goldstone mode, whose energy becomes 
nonzero only when lattice anisotropy is introduced. 

One might expect that the number of coupling constants might be so large that no useful 
information or test of the theory would be possible. As it happens, the fit to the energy 
of the gaps is overdetermined and the agreement between theory and experiment in some 
instances is quite remarkable, as can be seen in I. The observation of the modes whose energy 
depends on the in-plane anisotropy leads to the determination of the in-plane anisotropy of the 
exchange interactions. These quantities are difficult to obtain experimentally. Their values can 
be compared to calculationsEJll^ based on the electronic structure of the cuprates the knowledge 
of which may lead to a better understanding of the high-T^ superconductors. 

One should recognize that at the moment inelastic neutron scattering does not easily detect 
modes in the sub meV range of energy. As a result neutron scattering experiments have not 
detected those in-plane modes whose energy depends only on the in-plane anisotropy. Recently, 
however, the modes in the sub meV range of energy have been observed by ESR experiments 
of the group at RIKEN.tJO The mere existence of these modes tends to confirm the spin-wave 
calculations. Moreover, the fact that they are found in the predicted range of energy strongly 
supports the theoretical calculations in this paper. 

Briefiy, this paper is organized as follows. In Sec. II the Hamiltonian with its vari- 
ous anisotropic exchange interactions is specified. In Sec. Ill we start by discussing briefiy 
the framework within which the calculations are to be done and we give the Dyson-Maleev 
transformationtZl to boson operators. In Sec. IV the isotropic exchange Hamiltonian is dis- 
cussed, first within harmonic theory and then including spin-wave interactions, which are 
essential to obtain a qualitatively correct spectrum. In Sec. V the various anisotropics are 
included in an effective quadratic spin-wave Hamiltonian. In Sec. VI we give explicit results 
for the spin-wave energies for the case when the transverse wave vector is zero and show the 



comparison of our calculations with the recent experiments of the MIT group. In Sec. VII 
intensities of modes are discussed, with numerical results given for zero wave vector relative to 
the Bragg peaks for Cul and Cull. Our conlusions are summarized in Sec. VIII. 

II. HAMILTONIAN 

The Hamiltonian that we intend to treat is written as 

n = ni + n2, (2) 

where Tii includes almost all the significant interactions, namely all the intraplanar interactions 
and the unfrustrated interactions between nearest neighbors in adjacent CuO planes and 7^2 
includes small residual anisotropic interplanar interactions involving Cull spins. Since this 
latter term is totally negligible except for extremely small wave vector and for the lowest energy 
mode, it is only necessary to include contributions from 7^2 evaluated at zero wave vector. 
Since the effects of 7^2 are only relevant to the extremely low frequency spectrum, we defer 
consideration of 7^2 until Sees. V.4 and V.5. 
Thus we write 7Yi in tensor notation as 

7ii = 2 2^ Sj J/Sj + 2^ SiJi-iiSj + 2^ Sj J//Sj + 2^ JsSi ■ S^^i^^ , (3) 

{i,jei) {ieijeii) {i,j&ii) *g/ * '^''^ 

where i G I{i G //) means that site i runs over Cul (Cull) sites and ( ) restricts the summation 
to nearest neighbors of the indicated type in the same Cu-0 plane. The only unfrustrated 
coupling between planes is that (J3) between Cul's directly above or below one another. We 
will allow the couplings J/, J/-//, and J// to be anisotropic, whereas for simplicity we take J3 
to be isotropic. Here and below we use a hybrid notation for site labels in which the label i-\-r 
indicates a site at position r with respect to site i. In 7^2 we include the interplanar Cul - Cull 
and CuII-CuII couplings whose isotropic parts are frustrated. 

We first discuss the principal axes of the exchange tensor J/ associated with a bond between 
nearest neighboring Cul spins in a CuO plane. This bond is invariant with respect to two mirror 
planes: one in the CuO plane and the other perpendicularly bisecting the Cul - Cul bond in 
question. Accordingly, the principal axes of the Cul - Cul exchange tensor between nearest 
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neighbors lie along the three crystal (1,0,0) directions, just as they would be in the absence 
of the Cull's. In that case, the exchange tensor will have different values corresponding to 
the directions i) along the bond in question, ii) perpendicular to the bond in question but in 
the CuO plane, and iii) along the crystal c direction. The principal axes of the other in-plane 
interactions are similarly fixed by symmetry.EjO Then the Hamiltonian Tii may be written as 
follows 

iei 5i V / 

+ XI XI Jl-IlSiSi+52,1 + Jl-Il\^i ■ '^2,l][Si+52,i • h,l] + JJ'-lA^i ■ e2,l][Si+52,i • 62,1] 
ie//52,A ' / 

+5 E E [jhSISI^,^ + j\j[Si ■ 52][Si+s, ■ h] + Jh[^^ ■ e2][S.+52 • e2]) + J3 E Si • S. 

i&II (52 V / iel 



1 - , (4) 



where 5i {52} labels the nearest neighbor vectors in the plane connecting adjacent Cul's (Cull's) 
and (5i,2 labels vectors in the CuO plane which give the displacements of nearest neighboring 
Cul's relative to a Cull, and hat indicates a unit vector. Also ei, 62,1, and 62 are unit vectors 
in the CuO plane which are perpendicular to, respectively, 5i, ^2,1, and 82. 

We separate the Hamiltonian Tii into an isotropic part, Tio, and an anisotropic perturbation, 
Ti' . For that purpose we write 

AJi = i(jJ + jf)-J|, ^J,2 = Win + Jin) - J'Lii . AJ2 = liJh + Jh) - Jii , (5) 

(5Ji = |(J} - J/) , 6Ju = UJIii-JIii)^ 6J2 = Ujh-Jh), (6) 

J = ^{Jj + Jj + J'l) , J12 = -^{Ji-u + Ji-ii + Ji-ii) ' J2 = 3(^11 + Jii + Jii)- (''') 

Thus the AJ's describe the out-of-plane anisotropy (i. e. the energy which gives rise to an 
easy-plane) which is responsible for the 5 meV anisotropy gap in the spin-wave spectra of 
cuprates which do not have Cull's. Similarly, the (5J's describe the in-plane anisotropy (i. e. 
the anisotropy within the easy plane) and they i) are responsible for the weak ferromagnetic 
momenta's induced in the Cull subsystem by the staggered moment in the Cul subsystem and 
ii) contribute to the macroscopic or phenomenological four-fold anisotropy constant i^4.EJoQ 
(We shall see later that 7-^2 also contributes to K4.) Note that 6J12 is what was called Jpd in 



Refs. P and El but differs by a factor of two from its definition in Refs. O and 113. The largest 



coupling is J {J2I J ~ J12/J ~ 0.1 and J3/J ~ 10 ^), while the relative anisotropies, A J/ J 
and 5J/J are at most 10-3.|iiji3i§ 

With these notations the isotropic Hamiltonian is 



"^0 — 2Z^Z^ "^^i ' ^»+<5i ~^ A^ A^ Jl2^i ■ Sj+52 -^ + 2 Z^ Z^ "^28, ■ Sj+^j + 2_^ J^Si ■ S 



iel <5i ie// <52,i 

and the anisotropic perturbation is 



iell (52 



iG/ 



2+ o ci ' 



^' - -2^^1 zJ ^i^i+Si ~ ^Jl-i zJ ^i^i+52,1 ~ 2^J'i 2I ^i^i+S2 
iel,5i i&II,&2,i i€lI,S2 

_L 1 A 7 \^ I ex qy \ qy qx \ l^ j ST^ I qx qy , qy qx \ 

+ 2^-'l Z^ K'i '^i+5+ + '^J'^i+5+ I ~ 2"'^1 2^ K'i '^i+<5_ + -^i '^i+<5_ I 
je/,(5+ \ / ie/,(5_ \ / 

+SJ12 22 ^i^i+S^ ~ ^i^i+S^ 1 + '^-^12 zJ ^i^i+Sy - S^Sij^Sy 1 

je//,(5a; \ / ieii,5y \ J 






[Si ■ 82] [Sj ■ 82] - [Si ■ 62] [Sj ■ 62] 
where we introduce the following sums over the (5's: 



(9) 



5x = ±|ax , 5y = ±^ay , 5+ = ±|a(x + y) , (5_ = ±^a{x — y) , (10) 

as shown in Fig. ^. In Eq. (^, J = J + |AJ and similarly for the other J's. Since the 
anisotropy in the J's is so small (at most of order 10~^), we henceforth drop the tildes. 

It is convenient to express the spin components in a coordinate system in which one axis 
(the C, axis) lies along the line of the staggered magnetization. Thus we introduce the axes ^ 
and T) which are obtained from x and y by a rotation about the z axis of 7r/4. Then 

5- = (5« - s^'^)/V2 , sy = (5« + s^)/V2 , (11) 

so that 

JG/,<5i iG-f, 152,1 ieII,S2 

+ 2^Ji 2^ I Si 5'j_,_5^ — Si S^^g_^ I + ^6Ji 2_^ I Si 5'j^5_ — Si 5'j_^^_ I 

-SJl2 22 ^i ^i+S^ + ^i^i+5^ + '^-^12 22 ^i ^i+Sy + SiS^^g^ 
i€lI,Sx \ / iell,5y \ / 

-SJ2 E E Us] + s^s^) +6J2Y: E (sfs] + s^s^) , (12) 

JGe 5x'-j=i+25x V / *Ge (5a;:j=i+2(5y \ / 



where, in the last hne, i G e indicates that the sum is taken over only half the Cull spins, i. e. 
those on the e sublattice (see Fig. |1|). 

III. BOSON HAMILTONIAN 

A. Overview of the Calculation 

Since the Cul - Cull interaction is frustrated, the Cul and Cull sublattices are decoupled 
within mean-field theory or within harmonic spin-wave theory at zero wave vector. In other 
words, to calculate the energy gaps at zero wave vector we will need to include fluctuations, 
as first indicated by Shender.El Here, in view of the myriad of terms in the Hamiltonian, we 
need to proceed in as systematic a way as possible. In the original work of ShenderB it was 
found that the effective coupling between sublattices, which depends on fluctuations beyond 
mean-field theory or beyond harmonic spin- wave theory involved energies of relative order 1/S 
with respect to energies encountered in mean-field theory. Accordingly, here we will calculate 
all relevant effects in the spin-wave spectrum due to anharmonic perturbations up to first 
order in 1/5. Therefore we analyze perturbative contributions at one-loop order. To be more 
specific, we will introduce the usual Dyson-Maleev boson represent at iontZl of spin operators, in 
terms of which anharmonic perturbations involving three (four) boson operators are of relative 
order 1/Vo (1/5"). This means that we treat four-operator perturbations within first-order 
perturbation theory and three-operator perturbations within second-order perturbation theory. 
In technical language, this would be done by keeping all such contributions to the wave vector 
and energy-dependent self-energy. Since we work to low order, a more naive approach (which is 
entirely equivalent to calculating the self-energy) is both convenient and easy to follow. In this 
naive approach one truncates all four operator terms by contracting out pairs of operators in all 
possible ways. This reproduces exactly the results of the one- loop diagrams obtained by treating 
the four operator vertices in first order perturbation theory. In addition, we would note that all 
nonHermitian terms at order 1/S do not contribute to first order energies. So, at order 1/S we 
simply discard nonHermitian terms. Since the three-operator terms are of interest in producing 
small gaps, we will follow a calculational method which is strictly correct only at zero wave 



vector. The fact that in our treatment the small perturbations have the wrong dependence on 
wave vector is irrelevant because their effect is only nonnegligible very near zero wave vector. To 
avoid the algebraic complexities due to the fact that the magnetic structure has six sublattices, 
we simply construct, by the methods mentioned above, the effective quadratic Hamiltonian 
which includes all the self-energy corrections at order 1/ S. As a check that our calculations 
are really as consistent as we claim, we verify that the gaps have the expected dependence on 
the perturbations. In other words, when the perturbations are known to not produce gaps, 
our calculations reproduce that result. This type of check indicates that, for instance, our 
treatment of three-operator terms in second order perturbation theory is consistent with our 
treatment of four-operator terms in first order perturbation theory. 

B. Transformation to Bosons 

We make the following Dyson-Maleev transformationO to bosons (a, 6, .../): 



5+ = ^2Sa , 


S-- 


= V25aV(a) , 


Sl = S — a) a 


5+ = ^2Sb^ , 


s^ 


= y/2S(f){b)b , 


sl = -S + b^b 


5+ = ^2Sc^ , 


s; 


= V2^0(c)c , 


Si = -S + c^c 


St = V2Sd , 


sj- 


= V2Sd^(l)id) , 


sl = s - Sd 


S+ = V2Se^ , 


s: 


= V2^0(e)e , 


Sl = -S + e^e 


SI = V2Sf , 


sj 


= V2Sf<P{f) , 


S] = S-ff 



(13) 

where S^ = 5"'' ±iS^, 0(x) = 1 — x'^x/{2S), and we have left the site labels imphcit. In bosonic 
variables the isotropic interaction between spins assumes the form 

Saj ■ Sbj = si a\ai + bjbj + aibj + a -6] ^2 bjbjbjai + 6]a|ajaj + 2a]aibjbj 

Oai ' '^ej ^ J \ d^di ~r (^j(^j ~r diGj -\- (l^Cj \ 2 I GjCjCjdi -r CjCLj^d^di -\- Zdj^diCjCj 1 
Sai • S/j = Sl —didi — fj fj + difj + difj 1 ~ 2 [ '^ifj fj /?■ + fj'^i^i^i ~ "^(^i^ifj fj I 

Sbi ■ Sej = S I -bjbi - e]ej + b\ej + bic] \ - \\ bU]ejej + e]b\bibi - 26|6ie]ej j 
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S,, ■ S;, = Siblh + f]f, + h\f] + hfA - \ [h\f]f]f, + f0A + 26l6./j/,j 

The other interactions can be obtained by appropriate relabehng of boson variables. 
The effective bihnear spin-wave Hamiltonian is of the form (see below) 



^ = E 



A{c^),M^)^u^) + |i?(q)M4(q)^i(-q) + \B{c^)iM^)U-^) 



where ^i(q) = a(q) and so forth (in order 6, c, (i, e, and /). Here 

1 



«'■•'' - VN 



Ee"'"fi(q) 



uc q 



where N,,c is the number of unit cells. 



(14) 



(15) 



(16) 



C. Spin- wave Spectrum. General Considerations 

The transformation to normal mode operators Tfc(q) is 

J 3 

6(-q) = E%(q)M(q) + E^^.(q)^(-q) ■ 

3 3 

To preserve the commutation relations we require that 

P(q)pt(q) - Q(q)Qt(q) = J , P(q)QT(q) - Q(q)Pnq) = 

where X is the unit matrix. 

The transformation inverse to Eq. ([1^) is therefore 



r](q) = ^P,^.(q)et(q) -5:g,,.(q)a(-q) 

k k 

rM) = -EQ^.(q)*d(-q) + E^fc.(q)*a.(q) 



(17) 



:i8) 



(19) 



The equation that determines the normal modes is 



h(q),^]- = ^i(q)T-j(q) 



(20) 
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which gives 














A(q) B(q) 
_-B(q) -A(q)_ 




P.(q) 

_Q,(q)_ 


= ^j(q) 


"p.(q)" 

_Q,(q)_ 


where Pj is the column vector with components 


Plj,P2j,. ■ 


. Pnj and 



P = [Pi, P2, . . . Pn] and similarly for the Q's. 

From now on the arguments are always q. Then 

[A + B][P, + Q,]=cu,[P,-Q,-] 
[A-B][P,-Q,]=a;,[P, + Q,] 

Therefore 



f2r 



(22) 



[A + B][A - B][P, - Q,] = cu|[P, - Q,] . 



(23) 



Hence, the squares of the spin-wave energies are the eigenvalues of the matrix 



D(q) = [A(q) + B(q)]x[A(q)-B(q)]. 



(24) 



Roughly speaking the matrices A + B and A — B reproduce the stiffnesses in the two directions 
transverse to the sublattice magnetization. 

As we shall see later, for the Hamiltonian of the form oiTii these dynamical matrices assume 
the form 



A(q) 



an 


ai2C+ 


ai2C_ 





0156^ 


aieel 


ai2C+ 


an 





ai2C_ 


oiee; 


aisCy 


ai2C- 





an 


ai2C+ 


aiQCy 


ai^e; 





ai2C„ 


ai2C+ 


an 


aise: 


ai&dx 


aise* 


awCy 


aiee* 


ai^ex 


055 


„ Cx+Cy 

^56 2 


OieCa; 


ai5e*y 


OlsCj, 


oiee* 


„ Cx+Cy 

«56 2 


055 



(25a) 



and 
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B(q) = 



bii 


612C+ + 2J3SC, 


bl2C- 





bisSx 


foiee; 


612C+ + 2J3SC, 


foil 





fol2C- 


foiee; 


folSGy 


&12C- 





foil 


bi2C+ + 2 JsS'c^ 


biaCy 


foise; 





bl2C- 


612C+ + 2J3SCZ 


foil 


foise* 


foiee^: 


&i5e* 


bieCy 


foiee* 


foise^ 


fo55 


056 2 


bieSx 


bi5e; 


foisey 


foiee* 


fo56^ 


fo55 



(25b) 



where 



c+ = cos[a(g^ + gj^)/2] , c 



exp{iqya/2) , c^; = cos(g^a), Cj^ = cos(gya) 
cos[a(g^ - qy)/2] , c^ = cos(g^c/2) . 



(26) 



From now on we will analyze the energies of the modes for wave vectors of the form G + QzZ, 
where G is a reciprocal lattice vector. In that case the matrices A and B can be brought into 
block diagonal form consisting of three 2x2 blocks. The unitary transformation such that 
U'^^AU and U^BU are block diagonal depends on G, although, of course, the mode energies 
do not. For G = we have 



^ I/V2 1/2 
I/V2 1/2 



U 



1/2 
-1/2 

-1/v^ 1/2 -1/2 
-I/V2 1/2 1/2 
1/v^ 

1/v^ 



\ 



(27) 








1/v^ 
V 00 1/V2 -I/V2J 

U(G) for general G is given in Appendix^. For G = the transformed block-diagonal matrices 
corresponding to columns 1 and 2, (labeled "12"), those for columns 3 and 4 (labeled a = +1), 
and those for columns 5 and 6 (labeled a = —1) are 



A 



12 



A„ 



flu 
flu 

an + 2aai2 V^{ai5 + aaie) 
y V2(ai5 + aaie) (255 + o^ose 



(28a) 
(28b) 
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Bi2 = I , (29a) 

bn + 2abi2 + 2aJ^Scz V^ipi^ + afoie) , 

I • (29b) 

\ V2(6i5 + O-feie) &55 + 0-&56 



B. 



These results remain valid when 7^2 is included, providing it is evaluated at zero wave vector, 
which, as we have said, is an excellent approximation. 

D. Isotropic Interactions 

For a qualitative understanding of the mode structure we start by considering the results 
of linearized spin-wave theory when all exchange interactions are isotropic. Then one has 

an = 475" + 2Jj,S , aie = 615 = ^125* , 055 = 656 = 4725" , 612 = 275" (30) 

and all the other matrix elements are zero. 

In the "12" sector, we find two optical modes which are degenerate for all Qz, with 

{io/Sf = (4J + 2J3)2 - (2J3C,)2 ^ 16J2 . (31) 

Spin-wave interactions and anisotropic exchange interactions will have only negligible effects 
on these optical modes and accordingly we will generally not discuss these modes any further .Ej 
In the 0" = +1 sector we find modes with energies 

(iuX/Sy = 2J3(1 - c,) [8 J + 2 J3(l + c,)] ^ 16 JJ3(1 - c.) , 

iuf = . (32) 

Finally, the a = —1 sector has modes whose energies are 

{uj^/Sf = 2Jsil - c) [8 J + 2 J3(l + c,)] ^ 16 JJ3(1 - c) , 

c.^<^ = . (33) 

Note that all modes are gapless at zero wave vector and that for both a = +1 and cr = — 1 we 
have a dispersionless zero frequency mode due to the frustration of the Cul - Cull interaction. 
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Several aspects of the above results are noteworthy. First of all, as we will see from our 
calculation of the dynamic structure factor in Sec. VII, the a = +1 {a = —1) sector corresponds 
to modes in which the spins move out of (within) the basal plane and therefore we will refer 
to these modes as out-of-plane (in-plane) modes. (This identification can also be deduced from 
the way the mode energies depend on the out-of-plane and in-plane anisotropics.) For both out- 
of-plane and in-plane modes note the existence of a completely gapless mode: when the Cul's 
rotate in phase, they produce zero coupling on the Cull's, each plane of which can be rotated 
with zero cost in energy. The higher-energy out-of-plane and in-plane modes are degenerate 
because we have not yet included any anisotropy and these modes give rise to the usual twofold 
degenerate mode of the Cul subsystem. Even when more general anisotropic interactions are 
included, the higher-energy modes remain mostly on the Cul's and the lower-energy modes 
remain mostly on the Cull's. 

E. Mode Energies for General Interactions 

Here we give the mode energies in terms of the matrix elements of Eq. (|25|) for general 
interactions for wave vectors of the form q = (0, 0, g^). (The eigenvalues, but not the matrices, 
are invariant under addition of a reciprocal lattice vector G to q.) To evaluate Eq. (p^ ) within 
the low-frequency sectors a = ±1, we record the form of the two by two blocks. Since we need 
both A + B and A — B, we write 



[A + 7?B]„ 



ail + 2(Tai2 + r]bii + la-qJ^Sc^ + 2ar]bi2 


\/2[ai5 + o-ai6 + vhb + f^vhe] 


V2[ai5 + craie + r/615 + o-f/^ie] 


055 + 0-056 + Vhb + o-r?656 



(34) 



In evaluating Eq. ( ^4]) it is useful to note that in the a = +1 sector the matrix element 
[All + Bii]_|_ ~ 8JS is by far the largest matrix element. Similarly in the cr = — 1 sector 

[All 



Bii]_ ~ 8JS is by far the largest matrix element. In either case, then, Eq. (^) gives 
the squares of the mode energies as the eigenvalues of a matrix (or its transpose) of the form 



u 


V 




u 


V 


V 


w 


V 


w 



(35) 



where a/Z/m dominates all other matrix elements. In that case the eigenvalues are 
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{uj>f = Uu + 2Vv + Wv^/u, {uj<f = {UW-V^){uw-v^)/{u>f. (36) 

Explicitly, within the sectors a = ±1, we have 

U„ = au + 26i2 + a{2ai2 + hi) + ^J^Sc, , 
K = V2[ai5 + bie + o-(ai6 + &15)] , 
Wa = 055 + 656 + a(a56 + 655) , 
Ua = an - 2612 - 2J3SC;, + cr(2ai2 - &11) , 
v^ = V2[ai5 - &16 + o-(ai6 - ^15)] , 
Wa = 055 - he + cr(a56 - ^55) ■ (37) 



Substituting these evaluations into Eq. (^) [or, if need be, exactly implementing Eq. 
gives the four low energy modes for wave vectors along the c direction. Obviously, since the 
mode energies are derived from a two by two dynamical matrix, we can easily obtain exact 
expressions for their energies. 

IV. NONLINEAR SPIN WAVES 

A. 1/S Corrections to J, J3, and J2 

When we include the effect of spin- wave interactions at order 1/S on the Cul-Cul inter- 
actions or on the CuII-CuII interactions, we expect to get a simple renormalization. For the 
exchange interactions between neighbors in the same CuO plane, this effect is well known. As 
explained above, we decouple the fourth order terms in Sai ■ Sbj as 



— ^[hjhh'^i + hjolalai + 2alaihjbj] -^ —{alai + aj6j)[a|aj + hjh + o-ih + olh] (38) 

and those in Sgj ■ Sfj as 

- IH/J/]/, + fMe,e, + 2e|e,/;/,] ^ -(e|e, + e,/,)[ele, + /]/, + ej, + ej/j] . (39) 

From this result we conclude that J and J2 should be replaced by Z^J and Z2J2, respectively, 
with Zc = I — {1/S){alai + aibj), where i and j are nearest neighboring sites on the a and b 
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sublattices, respectively, and Z2 = 1 — (1/S')(e|ej + Cifj) in a similar notation, so that Z2 ~ Z^ 
at zero temperature. Zc has been calculated more accurately than this. (In Ref. ^ the value 
Zc ~ 1.17 is given.) For J3 we note that Oj and bj refer to sites in different CuO planes, in 
which case {aibj) ~ 0. So we should replace J3 by Z3J3, where 

Zs = l-{1/S){ala,) , (40) 

so that Z3/2 is essentially the magnitude of the zero-point staggered spin in the presence of 
quantum fluctuations. (Thus Z3 ^ 0.6 is very different from Zc.) 

B. The Effect of Spin- Wave Interactions on J12 

Now we discuss the effect of spin-wave interactions on J12, i. e. we consider the Shen- 
der interaction.! Correctly to order 1/S we construct the effective quadratic Hamiltonian by 
contracting two operators in all possible ways. I. e. we replace two operators by the thermal 
expectation value (indicated by (...)) of their product. Applying this procedure to the relevant 
terms in Eq. (0) we obtain the effective interactions between a Cul spin i on sublattice a and 
nearest neighboring Cull spins as 

Vae/Ju = a\ai { S - (aje+) - {e]ej) J + e]ej i S - (cjai) - (a|ai) J 

+ QiCj I S - {e]ej) - {a\e]) \ + a\e] I S - {a\a,j) - {aiCj) | (41) 

Vaf/Ji2 = aja, (-S - {a\i)) + {f]f,)\ + f]f, (s - (a jj) + (aja,) j 

+ a\f^ (s - (aja,) + {aj])\ + f]a, (s - {f] f,) + {alf,)) . (42) 

Here to leading order in 1/5 it suffices to evaluate the various expectation values with respect 
to the original quadratic Hamiltonian. At quadratic order we have symmetry such that (aje]) = 
{aiCj), {e]ej) = {fjfj), etc. We define 

jSs/Ji2 = S + (aifj) - (ajai) 

j'Ss/Jr2 = S+{a,f])-{f]f,) 

J12SIJ12 = S - {ttiCj) - {a\ai) 

j[?S/Ju = S- (a,e,) - (e]e,) . (43) 
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Note that j\^ - jjo^ = j\}) - j|?. Then 



'12 ^12 — ^12 ^12 



Vac — '-^12 ^^i^i ~r ^^19 ^^j^j \ '-'12 ^O^i^j \ '-'12 ^^{^j 

Vaf = -J^Saja, - j\SSf]f, + Ji?5aJ/, + jg^S/ja, . (44) 

Since we only work to order l/S*, we keep only the Hermitian part of these perturbations: 

t'ae ^ '-'12 '-"3'jQ'j + tV]^2 ^^j^j ' '-'12 '-' I '^i'^j ~r '-^j'^i I 

Vaf = -Al^Sala, - j\SSf]i) + J^2^sUi) + /ja.^j , (45) 

where 

j'^' = \[js+j?i]. '/r = ^['/ff +'/;?]. (46) 

As it turns out, the energies of the modes we study depend only on the single parameter 

« = (4? - J?^)S = {Jf^ - J^^)S = - Ji2((a.e,) + (a,/])) . (47) 

Note that the parameter 5 in Ref. ^ is 5 = a/ S. We evaluate this parameter in Appendix |B| 
and find 

a = C^JlJJ , (48) 

where C^ is a numerical factor which we found to be 0.1686. The anharmonic effects of Eq. 



(|45| ) give rise to contributions to the dynamical matrix of 

5aii = a , 

6aiQ = J 12 S — J12S , 

(5a55 = 2a , 

56l5 = '/|f5-Ji25. (49) 

It is knownl§ElB3 that in simpler problems these anharmonic effects give rise at zero momentum 
to effective biquadratic exchange interactions between sublattices which otherwise are frustrated 
in harmonic theory. To emphasize this point we treat a biquadratic interaction between nearest 
Cul - Cull neighbors (in the plane) which is of the form 
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^BQ = -^ E E(s. ■ s.+fe,j^ . (50) 

Then the contributions to the dynamical matrix are 

5aii = 4:JbqS , 

5ai6 = -2jbq5' , 

(^ctss = 8Jbq5' , 

6bi5 = 2jbqS . (51) 

Then using Eqs. (pTj) and (|36|) we find the mode energies at zero transverse wave vector (for 
large J) are now 

{uj>Y = 8JS[a + AjbqS + 2J;S{1 - c,)] = 8JS[a,s + 2J3^(1 - c,)], 

^2_4:aesJ3S{l-Cz){8J2S + aes) , . 

«eflF + 2J3b{l - Cz) 

where 

aes = a + 4:JbqS . (53) 

These results demonstrate that the Shender interaction does mimic a biquadratic exchange 
interaction at long wave length. However, in view of the relation for spin 1/2 that (Sj • S^)^ = 
Yq — |Sj ■ Sj, a biquadratic exchange interaction between two spins 1/2 is equivalent to a 
Heisenberg exchange interaction, and we may therefore assume that jbq vanishes. 

As before, there is degeneracy between in-plane and out-of-plane energies because we have 
not yet included anisotropy. However, by taking into account spin-wave interactions we now 
have the mode structure one would expect for an isotropic antiferromagnet: We have a doubly 
degenerate zero energy Goldstone mode at zero wave vector, and doubly degenerate nonzero 
energy modes for zero wave vector as shown in the right-hand panel of Fig. ^. The quantum gap 
in the optical mode uj^ at zero wave vector has been obtained for a number of other frustrated 
systems in several theoretical studies^'EiH^ beginning with the work of Shender.S However, 
because we have two subsystems which order at different temperatures, the emergence of this 
gap has a very unique signature not present in other experimental systems studied up to now.a 
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V. INCLUSION OF ANISOTROPIES 



A. Out of Plane Exchange Anisotropy 

To obtain the correct energy gaps at zero wave vector we must add the anisotropy due to 
anisotropic exchange interactions. (Since we are deahng with spin 1/2's, there can be no single 
ion anisotropy.) In this subsection we include out-of-plane exchange anisotropy. This part of 
the anisotropic exchange energy between sublattices a and b of the Cul's is given as 



Vab = -AJi 22 ^li^'bj^ij 5 
idajab 



(54) 



where Aj^ is defined so as to implement the nearest neighbor restriction. Thus, neglecting 
anharmonicity, we write 

Vab = 4AJ1 22 i^ai - Sai][Sbj " S^jj^ij 
ida^jab 

= IA.J,S J2 ia,-al){b]-b,)A,, 

i£a,jeb 



AJi^^[at(q)6(q) + 6t(q)a(q) - a\qp {-q) - a(q)6(-q)]^ 



(55) 



This result allows us to identify the contribution to the parameters of the dynamical matrix 
introduced in Eq. (|25|) as 



6ai2 = AJiS , Sbi2 = —AJiS 



(56) 



without having to explicitly consider the other Cul-Cul interactions. 

Next we consider the out-of-plane anisotropy of the Cul - Cull interactions. From the 
form of Eq. (^) we see that we only need construct the a-e and a-f interactions. For the a-e 
interaction we have 



K, 



- A J12 22 ^ai^ej^ij - 4 A J12 X] ( 'S'<i - 5'^i 1 ( S^j - S^j j A 

aGi,j£e \ 

|AJi25E 



ai -aUi e] - Cj A^- 



l^JuSY. 



q L 



a(q)e^(q)e^'i-("'=-"») - a^(-q)et(q)e^''-('"'=-'"'') + h.c. 
«(q)e^(q)e* - a^(-q)et(q)e*. + a\q)e{q)e^ - a(-q)e(q)e^ 



(57) 
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which gives a contribution to the dynamical matrix with 
Similarly 



Vaf='^AJuSY: 



a(q)/(-q)ex - a\q)f{q)el + a\q)f{-q)el - a(q)/^(q)e^ 



q 

from which we deduce that 



(59) 



(5ai6 = -|AJi2^ , She = |AJi2^ • (60) 

Finally we include the out-of-plane anisotropy of the CuII-CuII interactions. Thus 

v,f = -Aj2 E s:,s}^^., 

= 4AJ2 2_^ [S^i — S^i][Sfj — Sj-j]Aij 

i€e,jef 

= \AJ,S Y. (ej-e,) (/,-/]) A,, 

= AJ2SY.{e\q)f{q) - e\q)f\-q) - e(q)/(-q) + f\q)e{q)\ [c. + Cy] , (61) 

which leads to 

Sa^Q = 2AJ2S , 5b5Q = -2AJ2S . (62) 

The renormalization (at order l/S) of the out-of-plane anisotropy is accomplished by re- 
placing v/JAJi by Zg^/JATM 

It is instructive to see the influence of this anisotropy on the gaps at zero wave vector. 
Referring to Eq. ( PB| ) we see that the high energy mode gap due to the Shender fluctuation 
term, causes Uu to be nonzero. To check for gaps in the mode energies u'^ at zero wave vector 
it suffices to consider the quantity 

A = uw-v'^ = [2AJiS{l + a) + a] [2a + 2 A J25(l + a)] - [-V2aaf . (63) 

When we turn off both out-of-plane anisotropics, AJi and AJ2, the two modes u^ are gapless. 
When we allow the out-of-plane anisotropy to be nonzero, we clearly introduce a gap (A is 
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nonzero) in the out-of-plane (a = 1) sector but not in the in-plane {c = —1) sector. This result 
follows from the fact that the spins can still undergo a global rotation within the easy plane at 
no cost in energy. Hence we still have a single Goldstone mode with zero energy at zero wave 
vector. In order for this mode to have a gap, we have to take account of effects which lead to 
a fourfold in-plane anisotropy which we consider in the next subsection. 

B. In — Plane Exchange Anisotropy 

1. Cul-Cul Interactions 

In this subsection we discuss the effects of the in-plane anisotropy of the CuI-CuI exchange 
interactions. First of all, note that this perturbation is extremely weak. It gives rise to an 
effective fourfold anisotropy. This very small fourfold anisotropy only has a nonnegligible effect 
within the low frequency sector and even there only at zero wave vector. The Hamiltonian 
describing the in-plane anisotropy of the Cul-Cul interactions is 






(64) 



where j = i + 6, 6 is summed over four values (the two 5+'s and the two 5_'s), and ^(^-j-) = ±1. 
Then 



V,n = SJl E ^(^) 



i(ia,d;5 



(S - ala,)i-S + (3](3,) - i(25) U + «l0(a.) j ((3] + 0(A)/5, 



= SJi E (r{6)l-alai(3pj - lS[a, + aj][/?j + (3,] + \a\a\a0^ + /5,) 

i£a,d;5 \ 

+i(al + a^)p}pjpj - ^alaja^pjpjpA , (65) 

where ctj = a if site i is an a site and ai = d ii i is a d site, and similarly for (3j. We write 

Vin = ^"2,111 + ^"4,111 + Ki,!!! , (66) 

where the subscript 2 (4 or 6) indicates terms quadratic (fourth or sixth) order in boson oper- 
ators. Since we work systematically to first order in 1/S, we neglect V6,in- Also 
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iS 

= -6JiSj2 f [a^(k) + a(-k)] [6(k) + 6t(_k)] + [(;t(k) + rf(-k)] [c(k) + c^l-k)]") c+ 

<5,k \ / 

+5Ji55]('[at(k) + a(-k)][c(k) +ct(-k)] + [rft(k) + rf(-k)][6(k) + 6t(_k)]V_ (67) 

<5,k V / 



and 






(6J 



We now consider the effect of V2,in on the spectrum for k^ = ky = 0, so that c^ = c^ = 1. 
In this case because the perturbation is proportional to 6 — c or to 6"'' — c'l", one sees that V2,in 
only couples to the optical mode sector. Accordingly, we do not consider V2,in any further. 

We expect that this in-plane anisotropy should give rise to a macroscopic four-fold 
anisotropy. In order to obtain this anisotropy we must include anharmonic effects at rela- 
tive order 1/S. Now we decouple the four operator terms into quadratic terms times averages 
of the remaining quadratic factors. This calculation is done in Appendix y. In that calculation 
we naturally drop all contributions to the optical mode sector and of the rest keep only terms 
which have an effect on the mode energies at zero wave vector. The result is that contributions 
to the dynamical matrices due to the in-plane Cul - Cul interactions yield 

5au = IGC2T (69a) 

Sau = -4(6^2 - C2c - AC2b)T (69b) 

5hii = 8C2cr (69c) 

6bi2 = -lQC2bT , (69d) 



where r = {5JiY / J and the C"s are lattice sums defined in Eq. ( |C19|) of Appendix 0. It turns 



out that because r is so small, the only evaluation we need is that C2 = 0.01. Note that the 
contributions in Eq. (|69|) are of relative order 1/S which is consistent with the fact that they 
represent the effect of quantum fluctuations. The fact that they represent a modification in the 
zero-point energy is reflected by the appearance of the factor C2 ^ 1. 
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2. Cul- Cull Interactions 



Next we deal with the in-plane anisotropy of the CuI-CuII interactions. The terms in Eq, 
([T^) involving 5Ji2 are 



K5J12 = -5Jl2 Y^ 



■^i '~^i+5x ' i i+Sx 



ieii,Sy 






+Sy 



(70) 



In terms of boson operators this is 



Vsj,,=6JuJS/2j2 



[b — CjCiJ I Cli-\-x + C!'j_|_a.0(aj+x) + "i-x + di-xH^id-i-x) 



+ ( ej + 0(ei)ej j ( -25* + al^^ai+x + d\_^di_x j 



6J^JS/2J2 



[S - e\ei] h\ + <p{hi^y)hi-y + c|+ + 0(ci+j^)ci+j^ 



+ ( eJ + 0(ej)ei J ( 25' - h\_yhi_y - c{ 



+yCi+y 



+6JuJS/2j2 



[-S + fifi] tti-x + al_^(f){ai.x) + di+x + dl^^(j){di+x) 



+ U + fU{f^)] (-25 + aUa,_x + dlA+x 



-5J12V5/2E 



[S + flfi\ h\^ + (l){hi+y)hi+y + cj_„ + (t){ci.y)c 



-yj^t-y 



+ {h + fUih) 25 - blyh + y - CJ_„A 



y^l—y 



(71) 



This perturbation contains terms linear and terms cubic in the boson operators. The linear 
terms and (at relative order 1/S) the cubic terms will shift the equilibrium so that the boson 
operators are modified as 

Ci ^ ei + s , fi ^ fi + s , tti ^ tti + t , bi ^bi + t , Ci ^ Ci + t , di ^ di + t . (72) 

These shifts are evaluated in Appendix 0, where we find that (to leading order in 1/S) 
46Ji2\fs/^ . 2Jus J12SJ12JS/2 



t 



(73) 



SJ2 ' " 8J + AJ3 J2i8J + AJs) ■ 

These are the expected results. As one sees from Eq. ([12|) , the perpendicular field acting on 
an e spin is A6J12S in the positive rj direction, so that the perpendicular moment of the e spin 
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is ASe = 4:6Ji2Sxii = 4(5 7125*/ (8 J2), which agrees with v25's when Eq. ([73| ) is used. Further, 
due to the isotropic exchange, the field acting on an a spin is 2Ji2ASe = J12SJ12S/J2 in the 
negative rj direction. Thus ASa = —[Ji2SJi2S/J2]xi = —[Ji2SJi2S/J2]/[8J + 4:J3], which agrees 
with ^/2St when Eq. ([73|) is used. Note that A^e and ASa ai's both of order S, a result which 
indicates that the effects here are completely classical. 

To determine the effect of V^j^j o^ ^^e spin-wave spectrum we need to construct the effective 
quadratic Hamiltonian, which results from introducing shifts into anharmonic terms. This is 
done in Appendix 0. When we insert these shifts into the cubic terms of Vsj^^ ^^ ignore t in 
comparison to s because J ^ Ju- Thereby we get contributions to the dynamical matrix of 

Sa55 = 6aii = 25655 = 5JI2S/J2 = (S , 

5aie = -5ai5 = Sbie = -5bi5 = ^CS ■ (74) 

We also insert these shifts into the four operator terms of the isotropic Hamiltonian. As before 
we only keep terms arising from replacing two Cull operators by (e). The magnitude of other 
terms, e.g. Cul - Cul quartic terms when Cul shifts (a) are kept, are shown in Appendix 0to 
be much smaller than those we have kept. The result of the calculation in Appendix is that 
we get the contributions to the dynamical matrix of 

^055 = -(5" , 6b55 = -j(S , 

5a56 = -ICS , 6b,e = -CS . (75) 

Note that these perturbative contributions from the CuI-CuII in-plane anisotropy, are propor- 
tional to S, unlike the case for the other in-plane anisotropics. This indicates that the effect 
of 5Ji2 (which we called Jpd previouslyou) , is a classical effect which already appeared within 
mean field theory.H'a The other in-plane anisotropics only have an effect when we consider fluc- 
tuations. However, since the effect of 6J12 is rather small, we do not consider the effects of 
fluctuation corrections to it. 

3. Cull- Cull Intraplanar Interactions 

Here we consider the in-plane anisotropy of the interactions between pairs of Cull spins in 
the same plane. Their interaction is 
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V = -5J2J2 



ige 



7 . I ^i ^j + bi Sj \ + 2_^ y^i bj + ^i bj I 
=i+2Sx \ ) 5y:i=i+25y \ J 



■Sx-j=i+25; 



IS 



-^MhT.\is-eU 



i,Sx 






y:j=i+25y 

+ 



ej + 0(ei)ej 
ej + 0(ei)ei 






jjjj 



i,Sy \ 

5J2/572 [E Ue^ifj + /]) - (el + e.)/j/,] - ^ (4e.(/, + /j) - (eJ + e,)fjA 
^J^M^Pi^) ([/(k) + /n-k)]et(q)e(q - k) - [e(k) + et(-k)]/t(q)/(q - k)) , (76) 



where 



p(k) = -[cos(afc.)-cos(ay 



(77) 



This Hamiltonian is treated in Appendix 0, where the additional contributions to the spin- 
wave matrices (at g^ = 0) are found to be 



6a,, = -16[6Ji/J2][2C2a + ^2^] = -16^[2C2a + C25] , 6a,e = 16e[2C2a - Q 
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f75 



and 



Sh 



'55 



-16eC2fe, (5656 = 48eC26 , 



(79) 



where C2a and 6*26 are lattice sums defined in Appendix y. 

It is interesting to note that apart from a minus sign, these results are exactly the same as 
in Yildirim et al.li^ This difference in sign is to be expected because the Cull's are oriented in 
a hard direction with respect to only CuII-CuII interactions. Consequently, this term tends to 
decrease the gap. 

4- Cull- Cull Interplanar Interactions 

Here we consider the effect of interactions between a pair of Cull spins in adjacent planes. 
The situation we consider is shown in the left panel of Fig. H, where one sees that the isotropic 
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component of the CuII-CuII interplanar interaction is frustrated. To describe the anisotropy of 
this interaction we introduce the principal axes (shown in the right panel of Fig. ^) as follows 

ni = {—X + y)/y2 . 712 = {x + y) cos if)/ y2 + z sin il) ^ na = (x + y) sin-?/'/v2 — icos^ . (80) 

The angle ip is not fixed by symmetry. We then write the anisotropic CuII-CuII interaction 
HH'^' between nearest-neighboring spins i and j in adjacent planes asE^ 



^r'' = E^4s.-nF][s,-4^')], (81) 



A;=l 



where n^, is the fcth principal axes for the pair ij which can by obtained from the right panel of 
Fig. ^ by a rotation of coordinates, if necessary, and Kk is the associated principal value of the 
exchange tensor. The contributions of this interaction to the dynamical matrix are evaluated 
for q^ = Qy = in Appendix as 

5a55 = 6ae6 = 4(7^1 - i^ac^ - K3S^)S + 2{Ki + K^c^ + K2S^)Sc, , (82a) 

Sa^e = 6ae5 = 2(^2 - i^3)(c' - s^)Sc, , (82b) 

5655 = She = 2{Ki - K^s" - K,c^)Sc, , (82c) 

She = Sb,, = 2{K2 + K,)Sc,, (82d) 

where c = cosip and s = sinip. As we will see later, this interaction can only contribute 
significantly to the lowest energy in-plane mode, where its effect is through the combination 

Sia55 - ase + &55 - he) = ^S{K, - K^c^ - K:,s^){l + c,) = AAKS{1 + c,) . (83) 

Note that AK = for isotropic exchange. 

A closely related interaction is the long-range dipolar interaction, whose contributions to 
the dynamical matrix are also evaluated in Appendix 0. This interaction is dominant in 
Sr2Cu02Cl2". To include dipolar interactions we obtain (in Appendix ^) the result 

<5(a55 - 056 + &55 - ^56) = 65fVI'S'(l + Cz)X , (84) 

where X is the lattice sum 

X= J2 ^^^^ ' (85) 

jeir.Zij=c/2 ^ij 
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where i labels a fixed Cull site, aj is +1 if spins i and j are parallel and is —1 if they are 
antiparallel. Numerical evaluation yields 

X = 7x IQ-^A-^ . (86) 

Therefore we should replace AK by 

AKeff = Air + |//i|X. (87) 

5. Cul- Cull Interplanar Interactions 

Here we briefly summarize the results for a similar treatment of the CuI-CuII anisotropic 
interactions. The situation we consider is shown in the left panel of Fig. ^ where one sees 
that the isotropic component of the CuI-CuII interplanar interaction is frustrated. To describe 
the anisotropy of this interaction we introduce the principal axes for the CuI-CuII pair a — e, 
shown in the right panel of Fig. ^, as follows 

rhi = —y . 1712 = z cos (f) — x sin (f) , m-^ = —z sin </> — x cos <^ . 

The angle is not fixed by symmetry. We then write the anisotropic CuI-CuII interaction 
"HI^^^ between nearest-neighboring spins i and j in adjacent planes as 

^'' = E K[^i ■ ^SjW, ■ mSj)] , (89) 

fc=i 

where rhkiij) is the /cth principal axes for the pair ij which can by obtained from the right 

panel of Fig. ^ by a rotation of coordinates, if necessary, and K'^. is associated principal value of 

the exchange tensor. In Appendix we obtained the following contributions to the dynamical 

matrices for Qx = % = 

Sa,5 = 6he = l[K + ^2(1 - 3c') + K'^{1 - 3s^)] = Gj^n (90a) 

6aie = 6bi5 = j[k[ + ^2(1 + c') + K'^{1 + s^)] = Hlu , (90b) 

where c = cos and s = sin 0. We will see later that these terms have a negligible effect on the 
spin-wave spectrum. 
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VI. SPIN WAVE SPECTRUM 

Explicitly, the dynamical matrices corresponding to the effective quadratic Hamiltonian 
containing the above-mentioned anisotropics are of the form of Eq. (^) with 

an = AJS + 2 J3S + I6C2T + CS + a , 
ai2 = AJiS - 4(6C2 - C2c - 4C26)r , 

^16 = Jl2 S — 2 A J12S' + -^(S + Hj_ji , 

055 = 4J25 - 16^(2^2 - C2b) + 2a 

+4{Ki - K2C^ - K3S^)S + 2{Ki + K^c^ + K2s')Sc, , (91) 

a56 = 2AJ2S + 16^(2^2 - 3C26) - ICS + 2(7^2 - K^){^ - s^)Sc, , 

hi = 8C2cT , 

612 = 2JS - AJiS - 16C2bT , 

bi5 = J12 S — ^AJi2S — -^(S + Hj_ii , 

&16 = 2^'^12'S' + iCS + Gj^ll , 

&55 = ICS - 16^C2b + 2{K, - K2S' - K^^)Sc, , (92) 

he = 4J2^ - 2AJ25 -CS + 48^C2b + 2{K2 + K^)Sc, . (93) 

(In the above tabulation we not have included dipolar interactions. These are easiest to include 
when we give the mode energies because these terms can then be combined via Eq. (^) with 
the pseudodipolar terms which we treated explicitly.) 

In Table | we summarize the definitions of the various parameters and in Table ^ we give 
estimates of their numerical values. 

A. Cull's Ordered 

1. Without 1/5 Renormalizations 

Here we evaluate the energies of the four low-frequency modes in the presence of Cull 
ordering without any 1/S renormalizations. In what follows we will work to an accuracy of 
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about 1%. That is, the only corrections of relative order 1/ J we will keep are those of order 
Jvij J or J2I J ■ Then, in the notation of Eqs. p^ ) and ( |37D the components of the large matrix 

[A + aB]^ are 



V, = 8JS 



K = 2V2a JuS , 



W„ = 8J2S 



(94) 



We neglect terms which are small compared to a and obtain 



[A-Bl 



a=+l 



AAJiS + a + X3 



-V2> 



a 



-V2' 



a 



AAJ2S + 2a 



(95) 



for the out-of-plane sector, where X3 = 2J3S'(1 — c^), and 



[A + B],=_i 



CS + a + 6AC2T + X3 



V2{a - (S) 



V2{a - (S) 



2{CS + a)- 6A^C2 + AAK^sSil + c. 



(96) 



for the in-plane sector, where AK^s was defined by Eq. (|87|) . 
From Eq. (|36|) we get the higher frequency modes as 



(UJ 



[UJ 



>\2 



>\2 



8JS{AAJiS + a + zs) - SJia^a + 
8JS{a + X3) + 8Sa{-Ju + 2J2) - 



IGJoSa'^ 



a + AAJiS + X3 
16J2Sax3 



a + X3 



(97a) 
(97b) 



and the lower frequency modes as 



[UJ_ 



<\2 



6AJJ2S^i [4AJi5 + a + a;3][4AJ25 + 2a] - 2a'^]/{uj 



>\2 



(97c) 



uj<y = QAJJ2S' I [C^ + « + QAC2T + X3][2C^ + 2a- Q4^C2 + AAK^sSil + c,)] 

-2{a-csry{u:>r 

- / ^^'^■^^^'"^ [2x3 + 64(2r - OC2 + 8CS + AAK,sS{l + c,)] . 



(97d) 



In obtaining the above results we replaced UW — V"^ by UW with an error of order 1%. To 



obtain the last line of Eq. (|97d| ) we assumed that a dominates the other perturbations. 

As we have already seen, quantum fluctuations of the frustrated Cul - Cull interactions cause 
ui^ to be nonzero even if the exchange interactions are isotropic. When we introduce easy plane 
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anisotropy (by making A Ji and/or A J2 nonzero) we introduce a gap into cu^, but cjf has no gap 
yet, because without in-plane anisotropy a global rotation of spins within the easy plane costs 
no energy. The lowest mode develops a gap when we introduce the in-plane anisotropy and take 
account of quantum fluctuations. One might imagine that the strongest such anisotropy, namely 
that in J (scaled by the parameter 5Ji) would dominate in tuf . This effect is incorporated in 
the term proportional to r = 531/ J , and indeed when the Cull's are not ordered this term is 
the only one which contributes at g^ = 0. However, when the Cull's are ordered, the situation is 
different. Notice that this factor has no factor of S and more importantly, it is accompanied by 
the small numerical factor C2 ~ 0.01. These observations remind us that this effect is another 
fluctuation effect. Within harmonic theory or mean- field theory the anisotropy of these Cul- 
Cul in-plane interactions averages to zero. In contrast, the weaker in-plane interaction between 
Cul's and Cull's (scaled hy C, = {SJi2)'^/J) appears already in mean-field theory.a Thus, this 
term, which is proportional to S, has no factor analogous to C2 and it would dominate the 
term proportional to r except for the fact (see next section) that its renormalization factor 
Z(^ is quite small). However, when the CuH's are ordered, the interplanar CuII-CuII dipolar 
interactions contained in AK^s are dominant, and lead to the dramatic increase in the effective 
four-fold anisotropy observed at low temperatures. The isotropic interplanar nearest neighbor 
CuI-CuII are frustrated. The anisotropic CuI-CuII interlayer interactions (as embodied by the 
constants G and H) have only a negligible effect on the mode energies. 

2. 1/5 Renormalizations 

In this subsection we summarize how we incorporate the various renormalizations due to 
spin-wave interactions. We believe that the correct procedure is to calculate the mode energies 
correctly at first order in 1/ S and then set S = 1/2. Following this prescription we thereby 
obtain the following results 



(u;>)2 = 8JS 
(io>? = 8JS 



a + AAJ^SZl + X3^.? 



2 



^J^^Sa + ^^.yg^ „, (98a) 

a + 4AJi5z^ + 0:3^3 



a + x^Zl 



leJaSa^ 
a + X3ZI 



- 8 Jia^a + , ^ ^, (98b) 
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[UJ 



[UJ 



<\2 



<\2 



64772^2 ( [4A Ji^Z2 + a + x^Zl] [A^J^SZ^ + 2a] - 2J\ /{ujlf 

64 JJa^^ {[QS + a + UC2T + Zfxa] 

X [2C^ + 2« - 64eC2 + 4Ai^eff^(l + c.)] - 2[« - Csf\l{^^-? 
(UJJiS'^a 



(98c) 



[UJ 



>\2 



2Z^X3 + 64(2r - OC2 + 8CSZ^ + AAK^sS{l + c^)Z 



(98d) 



Here we noted that spins not in the same plane are essentially uncorrelated and hence we have 



Js ^ Z^Js , Ai^eff -- Z^AK, 



eflf 



(99) 



where Eq. ( 140|) gives Z3 ?^ 1 — 0.2/S' ^ 0.6. But since J3 and AKes always enter the 
mode energies in combination with an isotropic exchange constant, we associate with them the 
renormalizations 



-^3 ^ ^3-^3 , AKes -^ Z^AKcs , 



(100) 



where Z| = Z^Z^. Thus Z| = (1 - 0.2/5)(l + 0.085/5) = (1 - 0.115/5) -^ 0.77. Also, we 
will determine Z^ by comparison, in Eq. ( |108|) below, with the phenomenological treatments 
of the statics. For convenience we summarize in Table |T| the renormalizations of the various 
interactions which follow from our treatment to order 1/S. 



B. Cull's Disordered 



To get the energies of the spin-wave modes when the Cull's are disordered one sets J12 = 
J2 = (i. e. modes u^ and cuf no longer exist as elementary excitations) and a = 0, in which 
case we get 



[LUa 



8JS 



{uj_Y = 8JS 



AA.hSZl + 2.hSZl(l - c. 



64rC2 + 2J35Z|(1 



(101a) 
(101b) 



Note that in Eq. ( p7b|) we had dropped a term representing the four-fold anisotropy which 
is proportional to r, because such a term is negligible in comparison to a. Here, with a not 
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present, we restore this term in uj^. Note also that the higher energy mode is the one which 
has fluctuations out of the plane (as indicated by the dependence on AJi) and at zero wave 
vector is of the expected form, u"^ = 2HeHa, with the exchange field He = 4:JS and the 
anisotropy field Ha = 4AJiS'. The energy of this out-of-plane gap is about 5 meV in many 
lamellar copper oxide antiferromagnetsB The lower energy mode involves motion of the spins 
within the plane and would have no gap at zero wave vector except for the appearance of a 
small effective four-fold anisotropy, which was obtained previously^ from phenomenological 
considerations. The same result for the gap, namely u = 166Jiy/2C2S ~ 1.65Ji, is obtained 
from the microscopic calculation given in Appendix y and also in Ref. ^. 

C. Comparison of Static and Dynamic Theories 

Here we briefly compare our results with those of a mean-field treatment of the statics.Q 
In that calculation the four-fold anisotropy is included phenomenologically and the anisotropic 
Cul - Cull interactions are included even when the Cull sublattice is not antiferromagnetically 
ordered. When the Cull sublattice is ordered, the static treatment assumes that the Shender 
mechanism is strong enough that all spins are essentially collinear. So the dynamics of the Gold- 
stone mode should involve the static response coefficients, although spin-wave hydro dynamics^S 
rigorously applies only in the limit of zero frequency. 

Since the statics treat the four-fold anisotropy phenomenologically, as did Yildirim et al,c3 
we identify their four-fold anisotropy constant K, which scales the anisotropy energy per Cul 
spin, from 

E = -^K cos{Ae) , (102) 

because there are two Cul's per unit cell. Also 6 is the angle of the magnetic moment with 



respect to the easy, (1,0,0), axis. In Ref. |T3|the energy per Cul spin is (in the present notation) 

E = 22C2rS{SlSl/S^) . (103) 

So we make the identification K = 8C2TS, or, if we include the effects of the Cull's, 

K = 4C2(2r - ^)S . (104) 
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We start by comparing the results of the two approaches when the Cull's are disordered. 
There the spin-wave calculation completely ignores the presence of the Cull's, whereas in the 
statics the Cull's are characterized by their susceptibility in the pseudo-dipolar field caused by 
the small in-plane anisotropy of the Cul - Cull interactions. In the statics for temperatures 
far below the ordering temperature for the Cul sublattice (but still with the Cull's disordered) 
one has the effective fourth-order anisotropy constant fcgtat from the statics as 

fcstat = 2K + mpl^xAi - 8x7/^12]"' , (105) 

where we introduce the Cu spin susceptibilities, xi ~ 0.53/(8 J), xii ~ 0. 53/(8 J2), and (in the 
present notation) 

Mo = 4<5Ji2(S)x//, (106) 

where ( ) denotes a thermal average. If one takes 5Ji2 = 0.025 meV, then Mq = 2 x 10~^. 
Then the second term on the right-hand side of Eq. (|105|) is about 6 x 10^^ meV, compared 



with 2K which was founds to be 2 x 10~^ meV. So this correction (due to paramagnetic Cull's) 
which is absent from our spin-wave analysis is negligible. 

When the Cull's are well ordered, Ref. gives approximately 



''Stat 



2K + 8iSJi2f{Sf[0.53/{8J2)] . (107) 



Using Eq. ( p.04| ) as the identification of K, we see from Eq. ( p8cl] ) that the mode energy involves 
the combination (for ^ <^ r and (S <^ a) which we identify to be the effective value of k from 
the dynamics, kdyn, where 

kdyn = 8(2r - 0^2^ + CS^Z^ + AKesS^Zl = 2K + {5Ji2fS^ZjJ2 + AK^sS'^Zl . (108) 

We see that the term (0.53) (S*)^ in the statics appears as S'^Z^ in the spin-wave dynamics. With 
an appropriate renormalization Z^ ^ 0.19, these two terms are the same. Thus, as far as the 
intralayer interactions are concerned the comparison between statics and dynamics indicates 
that these terms are correctly treated. We also see that the treatment of the statics did not 
include the interplanar anisotropic interaction, A/Tefj. As we shall see, this term gives an 
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important contribution to the mode tuf , so it should be included in a reanalysis of the statics. 
In terms of the constant fcdyn we may write Eq. ( p8d| ) as 

(o;<)2(q = 0) = 64J2/.dyn ( j_j^^^^j ] ■ (109) 

Thus we conclude that except for the fact that the statics ignored the interplanar anisotropic 
CuII-CuII interactions, the two theoretical approaches are compatible with one another. In the 
next section will show that the experimental results from static and dynamic measurements are 
also consistent with one another. 

D. Comparison to Experiments 

The comparison between the present theory and experiments has been described briefly 
in several previous publications.Era'c3 Since a more detailed comparison is given in I, we will 
simply summarize the comparison of the theoretical and experimental results. First one has 
the estimate for J which is nearly the same for all cuprates. This estimate has been refined by 
Kim,c^ who gives J = 130 meV. The value J2 = 10.5 meV has been accurately determinedQ by 
comparing the experimental dispersion with respect to in-plane wave vector of Cull spin waves 
to various theoretical treatments which take account of spin-wave interactions.E°l 

Now we discuss the analysis of the magnon gaps at zero wave vector where values are listed 
in Table |IV|. We first fit the observedli^ in-plane gap when the Cull's are disordered. Equation 
([ToTbD yields uj_ = 16^/2C^6Jl and withlll u. = 0.066 meV, we get \6Ji\ = | Jy - J^|/2 = 0.042 



meV, a value which is about twice the theoretical estimates.o Using Eqs. (|104|) and (|105|) this 
corresponds to k = I6C2ST = 16(0.01) (0.5) (0.042)7130 = 1 x lO^^meV, compared to the 
value deduced from the statics,a k = 2 x 10^^ meV, for 70K< T < 120K. At low temperature 
(T = 1.4K) , where the Cull's are well ordered, the staticsEil gives k = 25 x 10~^ meV. 
From Eq. ( |109| ) with cuf = 0.15meV, we get A;dyn = 41 x lO^^meV. These results are listed 
in Table 0, where we see only a qualitative consistency between the interpretation of the 
static and dynamic experiments. It is possible that the quantum renormalizations (which 
affect the determination of k from the observed mode energy) are not quite correct. Also, the 
interpretation of the statics within which the CuII-CuII interplanar anisotropy is subsumed 

35 



into the four-fold anisotropy constant, k, is not strictly correct. If we fix 6Ji to fit the value of 
uj'2 at T = lOOK and assume that the interplanar CuII-CuII interactions result from the actual 
dipole-dipole interactions, then the temperature dependence of k results from the last term in 
Eq. (|108|) . With only dipolar (i. e. no pseudo dipolar) interactions, Eqs. (|87| ) and (|86D give 
(with g = 2.20) AK^s = 273 x lO^^meV, so that AK^sS^Zi = 53 x 10"^ meV, from which 
^dyn = 56 X 10^^ meV. From Table |V| it is clear that the experimentally deduced temperature 
dependence of k is qualitatively accounted for by the intraplanar dipolar interactions, especially 
if one increases AK by assuming it to have a small pseudodipolar component. 

Now we consider the higher energy modes. Fitting to the observed^ energy u;_|_ = 5.5 meV 
of the out-of-plane gap when the Cull's are disordered to Eq. ( |101a|) (with Zg = 0.6) we obtain 
the value of AJi = 0.081 meV. As was the case for 6Ji, this result is also about twice the 
theoretical estimates for a simple CuO plane.t^'tJ Given the values of these parameters, both 
higher-energy modes at low temperature involve only the one additional parameter, a. If we 
determine a from u^ we get a = 0.14 meV, whereas if we determine a from u^ we get a = 0.13 
meV. These two values agree perfectly with one another and their average coincides with the 
theoretical evaluation of Appendix that a = 0.13 meV. Clearly these agreements strongly 
support our interpretation of the role of fluctuations embodied by the parameter a. Note that 
a biquadratic interaction between two spin 1/2's can be subsumed into a ordinary Heisenberg 
exchange interaction. Therefore biquadratic exchange can not contribute to Ai^gff- 

Finally we consider the lower energy out-of-plane mode in the zero temperature limit. The 
AFMR datap3ll3 gives cu^ = 1.7473(4) meV, more accurate than, but entirely consistent with, 
the data of Ref. |^. Evaluating the expression in Eq. ( |98c|) with AJ2 = gives u^ = 1.717 
meV. If we fix AJ2 to fit the experimental value of this gap, we get AJ2 = 0.004 ± 0.004 meV. 
We attribute a large uncertainty to AJ2 because its value changes significantly if AJi or a is 
slightly modified. To get the same relative out-of-plane anisotropy, A J/ J, for the CuII-CuII 
exchange as for the Cul-Cul exchange would require AJi = 0.008 meV. 
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VII. DYNAMIC STRUCTURE FACTOR 



The cross section, cr(q, uj), for inelastic neutron scattering from magnetic ions is proportional 
to the dynamic structure factor ^"^(q, cj) which in turn is related to the spin-spin correlation 
function. We have 



a/3 



(110) 



According to the fluctuation-dissipation theorem, we may write 

S'"^{q,u}) = -n(cu)Imx"^(q,cj-iO+) , (111) 

n 

where n{uj) = [e^'^Z^'^^) — 1]^^ and, in the usual notation,E3 the A — B Green's function is defined 



as 



{{A;B))^ = Y.P. 



(?2|A|r7i)(r7i|B|n) {n\'B\'m){'m\A\n) 



^ — Em + En 



U + Em — E. 



m -^n 



;ii2) 



where \n) and \m) are exact eigenstates with respective energies En and Em and pn is the 
Boltzmann weight of the state \n). Then x, the dynamic susceptibility, is written as the 
Green's function 



X'^'^(q,c.) = ((S"(q);5^(-q)))_ 



;ii3) 



We construct the dynamic susceptibility by writing the spin operators in terms of boson oper- 
ators at leading order in 1/5: 



(114) 



5+(q) = V2S (a(q) + 6t(_q) + ct(-q) + d(q) + et(-q) + /(q) 
S'i-q) = V2S (at(q) + 6(-q) + c(-q) + rft(q) + e(-q) + /t(q) 



Thus we have 



5^(q) = [5+q) + 5-(q)]/2 = ^^J] [Vmiv)U^) + Kn(r/)*d(-q) 



(115) 



and 



s'ici) = -z[5+q) - s-{q)]/2 = ■s/s/2Y: [Vm{z)U^) + K.(^)*d(-q)] , (US) 
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where the operators are labelled as in Eq. ([T5|) and the transpose of the column vectors V(q;) 
is 



y(r/) = (l,l,l,l,l,l), 



V(z) = 2(1, -1,-1, 1,-1,1). 



;ii7) 



Thus we may write 



X 



a/3 



;q,u;) = iS5:(([Kn(«)e™(q) + Kn(a)*d(-q)] ; [K(/?)en(-q) + K(/5)*ei(q)])).. ■ (US) 



We may evaluate these response functions in terms of normal modes. Suppose we have 
found the unnormalized right eigenvectors of the dynamical matrix, Eq. ( ]23|) . That is we have 
the column vectors $j which satisfy 



[A + B][A-B]$,- = ^,^$ 



J - ^j ^j 



;ii9) 



Then we make the identification that 



Pi - Qi = ^j^j ■ 



(120) 



We can arbitrarily fix the phase of the normal mode operators so that Xj is real positive. Then 



[A - B]a;,<l>, = [A - B][P, - Q,] = ..,[P, + Q, 



■jj ' 



'1211 



or 



so that 



P. + Q. = (^.M)[A-B]<|., 



3 1 



p — ^ 
'~ 2 



O- = ^ 
^j 2 






$. 



To use Eq. (p^ ) we write 



PlP,-Q]Q, = |($l[A-B]$,) 



so that 



'122) 



:i23) 



(124) 
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Ui 



X, 



N[A-B]$,. 



(125) 



Then we write the susceptibihties as 



(2/5)x"''(q,cu)= E [K«(«)P™r(q) + Kn(a)*g™r(q)] 

•m,n,r 

X [VnWPnMT + K(/3)g„r(q)1 ((r.(q); rt(q))), 
+ Y^ [Kn(a)<5mr(q) + V'm(a)*Pm,.(q)] 

mnr 

X [K(/5)*gnr.(q)* + V;(/9)P™.(q)1 ((rt(q); r,(q))). 



E< 



VMP, 



V(a)*Q, 



V(/?)PJ+ V(/3)*Q, 



[cj-cj^(q)] 



-1 



+ 



V(a)Q, 



V(a)*P, 



V(/?)QJ+ V(/5rPr 



[cj + c^^(q)] 



-1 



E 



Jf{^) ^ If{ci) 



(126) 



where we left the argument q imphcit in several places. We will refer to I and J as 'intensities', 
although to get inelastic neutron scattering cross-sections one needs to include several other 
factors. At low temperature we only need 



/r (q) = a: 



X 



(5«,J\/(^)t$,(q)j+cu,(q)-X,JV(^)^[A-B]$,(q) 
(5^,, f V(z)t$,(q)^ + uj^yHp,^ {v{n)\A - B]$,(q) 



(127) 



In writing this result we used the fact that V(?7) is real and V(z) is imaginary. From now on, 
we specialize to the case of wave vectors of the form q = G + qzZ. In that case /^^ + J^'' vanishes 
and 



V(^)t<I.,(q) 



^r(qj 



:i28a) 



$,(q)t[A-B]<I.,,(q) 



V(r7)t[A-B]<I>,(q) 



rm 



(128b) 



cc;,(q) U,(q)t[A-B]<l>,(q) 
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The above results are useful for the out-of-plane (cr = +1) modes in which case [A — B] is the 
small matrix. Alternatively, for in-plane {a = —1) modes when A + B is the small matrix the 
following forms are useful: 

\ 2 



TZZ 

r 



V(z)t[A + B]^,(q) 



:i29a) 



u;,(q)U,(q)t[A + B]^,(q) 



ir 



V(r7)tvl/,(q) 



UJr{q) 



(129b) 



^,(q)t[A + B]^,(q) 



For high symmetry directions of the wave vector, the matrices A and B may be brought into 
block diagonal form by a unitary transformation U. In that case we may apply the above 
formulas in terms of the transformed quantities indicated by primes: 



A' = U^AU , B' = V^BV , 

$; = U^$, , ^; = U^^, , V'(a) = U^V(a) 



(130) 



For wave vectors which are equal modulo a reciprocal lattice vector, the corresponding 
quantities. A', B'. $', and \1'' are equal. However, the intensities at such equivalent points 
will differ because U, and hence V, depend specifically on the zone of the wave vector. This 
can be seen explicitly in Appendix ^ where we obtain the results summarized in Tables ^ 
and |VI1| . Note that the a = +1 sector does have intensity mainly in P^ in confirmation of our 
identification of this as the out-of-plane sector. Similarly, the a = —1 sector has its intensity 
mainly in I^"^ as expected for in-plane modes. These identifications are also consistent with 
the fact that the a = +1 modes depend on the out-of-plane anisotropies scaled by the AJ's, 
whereas the a = —1 modes do not involve these quantities. 



VIII. CONCLUSIONS 



Here we briefly summarize the significant conclusions from this work. 
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• 1. The degeneracy, present within mean- field theory, in which the Cull sublattice spins can 
be globally rotated with respect to the Cul spins is removed by quantum fiuctuations which 
cause the sublattice magnetizations to be collinear, as first indicated by Shender.H 

• 2. A degeneracy present within mean-field and linear spin-wave theories, in which the mag- 
netization can be globally rotated through an arbitrary angle within the easy plane is similarly 
removed by quantum fiuctuations, as first proposed in Ref. |1^. 

• 3. These fiuctuation effects, in addition to selecting the ground state from among the classi- 
cally degenerate configurations, also give rise to nonzero energies of the corresponding spin-wave 
excitations. The most dramatic evidences of this phenomenon are the striking increases of the 
out-of-plane gap energy from 5 to 10 meV and that of the in-plane gap from zero to 9 meV 
when the Cull sublattice evolves from disorder to order. 

• 4. The experimental results of inelastic neutron scattering for the lowest energy gaps are 
broadly consistent with the effective four-fold anisotropy previously obtained from the statics 
experiments.u More precise agreement may depend on more accurate understanding of the 
various renormalizations due to quantum and thermal fluctuations. 

• 5. Our improved theoretical treatment which now includes the interlayer dipolar interactions 
resolves the mystery surrounding the dramatic increase (flrst found in the staticso) in the 
effective four-fold anisotropy as the temperature is reduced into the regime where the Cull's 
order. In fact the dipolar interlayer interactions between the Cull's dominates the effective 
four-fold anisotropy when the Cull's develop long range order. 

• 6. Recent AFMR resultsi^ lead to an identiflcation of the small in-plane anisotropics and 
qualitatively conflrm previous theoretical estimates of the exchange anisotropy induced by spin- 
orbit interactions.EJo 
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APPENDIX A: INTENSITY CALCULATIONS 

In this Appendix we evaluate the intensities for which formulas are given in Sec. V. We first 
give the unitary transformation which brings the matrices A and B into block diagonal form. 
We do this for wave vectors q = G + g^i, where 



G = 27r 



Hx Ky Lz 
a a c 



(Al) 



where H and K are either both half integral or both integral and H + K + Lis an even integer. 
Then 



/ 



U = - 
2 



V 



V2 










^/2 
-1)2^ v^ 









1 













1 











\2H 



\ 



^ 



(A2) 



(2)2^ V^ -(2)2^ V^y 



The first two columns are the high frequency Cul optical modes. Columns #3 and 4 are the 
a = 1 out-of-plane modes and columns ^^5 and 6 are the a = —1 in-plane modes. The following 
results hold for all wavevectors of the form q = G + QzZ. 

1. Out-of-Plane Modes 



For the out-of-plane sector we have (for dominant J) 
3:3 + AAJiS + ift a/V^ 



A'-B' 



2~ "/ V - I ^ A' + B' 

a/V2 4:AJ2S + a 



8JS V2J12S 
y2Ji2S 8J2S 



(A3) 



independent of G, where X3 = 2J3S'[1 — cos(cg2/2)]. Note that q^ is measured relative to 
the reciprocal lattice vector in question. We now tabulate the right eigenvectors of the block 
matrices M_|__ = [A' + B'] [A' — B'] associated with the eigenvalues (the squares of the mode 
energies) u^. We have 
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$> = [1,0], {u 



>\2 



(8JS)(x3 + 4AJi^+ia 



2"y ) 



$< = [-a/V2, X3 + AAJiS + la] , (cu<)' = {8J2S) 
Also we find that 



AAJ2S + a- 



\a^ 



X3 + 4AJi5 + |a 



Vf^)' 



^iz)' 



"l-(-l 


^i^+i<' 






- 


z2^V2 





V(r^)' 



1 + [-l)H+K 



for integer H 



V(r^)' 



for half — integer H 



(A4) 



(A5a) 



(A5b) 



Note that the vectors V(a)' depend on G. Substituting these evaluations into Eq. ( |128| ) we 



obtain the results for the intensities in Tables ^ and [VI H for the out-of-plane (a = +1) modes. 



2. In-Plane Modes 



For the in-plane sector we have (for dominant J) 



A' + B' 



a;3 + C5+|« V2{CS-\a) 



A'-B' 



SJS V2Ji2S 
V2JuS 8J2S 



(A6) 



and we now tabulate the right eigenvectors of the block matrices M_+ = [A' — B'] [A' + B'] 
associated with the eigenvalues (the squares of the mode energies) u"^. For dominant J we have 
the approximate results 

[1,0], (cu>)2 = (8J5)(x3 + C^+i«)' 



^ 



> 



[UJ 



<\2 



-v^(C5-ia),X3 + C5+i 



a\ 



{SJ2S) 



2CS + a-2 



(xs + CS+^a) 



and 



y(z)' 



Viz)' 



-v^(-l)^ 

0^ 




V(r^)' 



1 _ (-l)H+K 





for integer H 



nvY 







-V2{t 



x2H 



for half — integer H . 



(A7) 



(A8a) 



(A8b) 



As before, only the vectors V(a)' depend on G. Substituting these evaluations into Eq. (|2 



we obtain the results for the intensities in Tables ^ and |VI1| for the in-plane (cr = —1) modes. 
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APPENDIX B: SHENDER PARAMETERS 



In this Appendix we evaluate the averages 

Ai = {amfl) 

M = {amCn) 



(Bla) 
(Bib) 



where site n is a nearest neighbor of site m. The above quantities can be calculated perturba- 
tively in the frustrated coupling J12 between Cul's and Cull's. (See Fig. |I].) 



1. Ai 



Thus 



A, = -(0 



Vl-II-ttrafl 



0-0 



amfljVi-n 







(B2) 



where S is the unperturbed energy of the virtual state relative to the ground state. Here we 
invoke perturbation theory relative to decoupled Cul and Cull subsystems, and 



Vi- 



II — -'12 



J12O 



J2 [a-hi + e]ej + a^ej + aje]] + J2 ["-hi + fjfj + 4/j + /] 



i(^a,5 



iGa,5 



+ Y: [blb^ + e]ej + h\e, + e]h] + ^ [h\h + /]/, + hf^ + h\f] 



i&b,S 



ieb,S 



+ Yl [4ci + e]ej + cjej + e]ci] + ^ [cJq + f]fj + Cifj + cj/ 



"ii 



i£c,(5 



j£c,(5 



+ Y. [dUi + e]ej + diCj + d\e]] + Y l^Ui + f]fj + 4/i + f]di. 
Only terms in Vj-n which have operators in both subsystems contribute. So, effectively 



(B3) 



i,5 



ttiCj + a|e] + alfj + /ja^ + bjej + e]bi + kfj + h\f] 



+c\ej + e]ci + Cifj + cj/j + djej + d\e] + 4/i + f]di 



(B4) 



where site j is the appropriate nearest neighbor of site i. In fact, in Eq. (p^ ) we need to have 
only terms with / or e^ and a\ d\ b, or c. So we set 
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v7-// = Ki = Ji2^E 



e] + ajfj + e]hi + hfj + e]ci + q/j- + d\e] + d\fj 



(B5) 



Thus with n = m + 5af we have 





-ji25>;(o 


-Ji25>;(o 


-Jl25)^(0 


-^12^ > ; (o 


-Ji25>;(o 


-^12^ > ; (o 


-J12S > ; (0 



1 



'^m j-7 Ctj 



1 



'£ 



-.4 



0)(0|/i,+,„,[eU, + /.+5„,]|0) 



(B6) 



Here we neglected the energy of the Cull modes in comparison to that of the Cul modes. Also 
we used the unusual notation that 



6st = rt-Ys . 



(B7) 



Then 



Ai 



J12S 






EE(o 



a'lqj^aiqj 



0) e**^''^''-'''"'' 



X 



gt(k)eik.6„e ^ f(-]^)e^^-^-f]p(-]^)e-'^-^-f\0) 



uc q,k ieb 



b{-^)ga{q] 



) e**^''^'^-'''"'' 



gt(k)e«k.5,, ^ /(_k)e*'^-^'^]/1'(-k)e-*'^-^"^|0) 



uc q,k iec 



ci-qj^a(qj 



q\ g*(q+k)-ri; 



q\ gi(q+k)-ri™ 
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X 



(0|[et(k)e*^-'''*- + /(-k)e*'^-''*]/t(-k)e~*'^-''-^|0) 



JuS 






EE 



aiqj^a'(qj 



q\ g«(q+k)-r,. 



x(0|/"f(-k)e-^''-^«^[e^(k)e^'^-^- + /(-k)e*^-'''"-f]|0) 



JloS 



EE 



a(q)^^(-q; 



0) e**^'^^''^''"*' 



Ar2 

uc q,k ieb 

x(0|/^(-k)e-*''-^«^[e^(k)e^'^-^'"= + /(-k)e*'^-''''^]|0) 

JloS 



EE 



Ar2 

uc q,k lec 



aiqj^c(-qj 



q\ g*(q+k)-ri; 



x(0|/^(-k)e~*'^-^«^[e^(k)e*'^-^- + /(-k)e*^-''-^]|0) 

q\ gi(q+k)-ri™ 



-%^i:i:(o«(q)^<''(q) 

^^uc q,k ied ^ *-' 

x(0|/t(-k)e"^''-^«^[e^(k)e^''-^* ^ f(-\,U^-^^f 



/(-k)e^ 



(B8) 



Doing the sum over i we get 



A, 



J12S 



N, 



J12S 



uc q \ 

E(o 



J12S 
aT 
^12 



Eo 



'uc q 
J,2S 



a iqj^aiqy 



K-q)^a(q; 

c(-q)^a(q^ 



.r E(o^^(q)7«(q^ 

uc q \ ^ 



Jl2^ 

iVuc 

J12S 
J12S 
h2^ 



so 

q ^ 

E(o 



E(o 



'uc q 
J19O 



«(q)^«Hq) 



a(q)^&(-q' 



aiqj^c(-qj 



0\ (0|[e^(-q)e~''i-''- + fiq)e-'''-^-f]f{q)e'''-^-f\0) 

o\ (0|[e^(-q)e-*'i-''^- + /(q)e-^'i-^''^]/t(q)e*'i-^"^|0) 

0) (0|[et(-q)e-''-''- +/(q)e-''-'=/]/t(q)e^'i-^"/]|0) 

o\ (0|[et(-q)e-^'i-^* + /(q)e-^'i-^*]/^(q)e^'i-''«/|0) 

0) (0|/^(q)e''i-''"^[e^(-q)e-'^-^- + f{q)e-'''-^-f]\0) 

o\ (0|/t(q)e*'i-'''"^[et(-q)e-*'i-''''^ + /(q)e-*'i-^*^]|0) 



AT,, 



^ 



E(0 «(q)7^^(q) 0) {0\f\q)e'^-'-f[e\-q)e-''^-'-^ + f{q)e-'^-'-f 



For the Cull subsystem we have the usual relations 



(B9) 



e(q) = M(q) -m^6^{-q) 



/n-q) 



-mqr7(q) + /q5^(-q) 



(BIO) 
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where ri{q) and 6{q) are the normal mode operators for the Cull subsystem and /q and rriq are 
given in Eq. ( pi)ll| ). In the ground state we evaluate the averages to get 



Ai 



J12S 

iVuc 

J12S 

J12S 
A^ 
J12S 
A^ 
J12S 

'a^ 

J12S 

A^ 
J12S 

"XT 
J12S 

A^ 
J12S 

J12S 
A^ 
J12S 

"aT 

h2 



>;(o 

q ^ 


>;(o 

q ^ 


>:{o 

a ^ 



aH 


q)^a 


(q) 


b{- 


.)i 


a(q) 


c(- 


-)i 


a(q) 



) [-/qmqe-^'i''- + /2e-^i-^«/]e^'i-^«/ 



) l-l^m^e-'"^-^'^' + /2e-iq-V]e^q-^-/ 



0) [-/qmqe-^''''- + lle-'''-^^f]e'''' 



E(o ^nq)7a(q) o\ Hqmqe-*'i-^^'= + /^e-*'i-^*]e''i''"^ 



>:(o 

q ^ 


>;(o 

q ^ 





E 



£ 



aiqj^a'(qj 



«(q)^K-q) 



alqj^c(-qj 



0) e*'i''"-f [-/qmqe^^'i-^'"^ + mje"*'^'^"^ 






E(o «(q)7^Hq) o\e*'i-^"^[-/qmqe-*'i''*+m2e-''i''''^] 



E(o 

q ^ 

E(o 



E(o 



>uc q 

J12S 



£ 



a"^(q)^a(q) 
^(-q)^a(q) 



c{,-(l)-a[q} 



0) [-/qmqe^'^'?^ + /^] 

0\ [-/qmqe-*'^"/2 + ;2gig,a/2|gig.a/2 



) [-/qmqe*«^"/2 + /2e-ig.a/2]gig.a/2 



^^^ E(0 d\q)-a{q) ) [-/qmq + /^e^"-] 



J12S 

A^ 

Jl2^ 



E 



AT, 



E(o 



uc q 



J12S 

'aTu 

^^12^ 



E(o 



•uc q 
J-\oS 



a(q)^a"^(q) 
a(q)^&(-q) 



alqj^c(-qj 



0)Hqmqe*"''^+m2] 



) [-/qmqe-*^""/' + m'e 



iqya/2 , 2 igya/2] iga;a/2 



o\ [-l^m^e"iy''l^ + m^ 



g Jgya/2jgii}a;a/2 



A^„ 



5:(0a(q)-dt(q 



£ 



) s - /qiTiq + m^e 



2Jqxa.] 



(Bii; 



For the Cul subsystem we have normal modes via the transformations, 



a(q) = (l/V2)[a+(q) + a_(q)] 
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rf(q) = (l/v/2)[a+(q)-a_(q)] 
%) = (l/v^)[6+(q) + 6_(q)] 
c(q) = (l/v^)[6+(q)-6_(q)]. (B12) 

In terms of these operators (in the order a+, 6+, a_, 6_) we have the matrices A and B: 



A(q) 
S 



AJ + 2 J3 














4J + 2 J3 














4J + 2 J3 














4J + 2 J3 _ 



and B(q)/S' as 






2J(c+ + c_) + 2J3C^ 












2J(c+ + c_) + 2J3C^ 


























2J(c+ - 


-C_) + 2J3C^ 








2J(c+ - 


-C^) + 2J3C, 








(B13) 



(B14) 



where 



c+ = cos[a(ga; + qy)/'A ' c_ = cos[a(g^ - qy)/2] , c^ = cos{q^c/2) 



(B15) 



Now each sector has relations analogous to the Cull's: 

aa(q) = ^a,qtta(q) - m^,q/5|(-q) 

where Q;o-(q) and /3o-(q) are the normal mode operators, and 



Here 



cr,q 



^ + EM) 
2i?.(q) 



m 



cr,q 



2i?.(q) 



^cr,q''^cr,q 



i?.(q) 
2i?.(q) 



(B16) 



(BIT) 



i?.(qf = A2-i?.(q)^ 



(B18) 



where 
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so that 



Thus 



A = 4J + 2 J3 
B^{q) = 2j( cos[{qx + qy)a/2] + o-cos[(g^ - qy)a/2] j + 2 J3 cos(g^c) , (B19) 



i?+(q) = 2J3COs{qzc) + AJ cos{qxa/2) cos{qya/2) (B20a) 

B_ (q) = 2 J3 cos(g^c) - 4 J sin(g^a/2) sm{qya/2) . (B20b) 



(0|at(q)ia(q)|0) = ^ E(0|4(q)«.(q)|0)£^.(q)-^ 



Similarly 



(0|a(q)|at(q)|0) = ^ E(0|<'"(<l)4(q)|0)-B,(q)-' 



±Ee,^,«l)-'-i:^J^ (B22) 



(0|<it(q)^a(q)|0) = ii:a(0|4(q)a,(q)|0)B,(q)-' 



25?-"-^'«')"-?^15W '"''' 



{0|a(q)irf'(q)|0) = ^ E<'(0|«-(q)«i(q)|0)K(q)-' 



^E^^.W-'-E-^^ (B24) 



(0|i(-q)ia(q)|0> = ^ i:(0|6„(-q)a„(q)|0>B„(q)-' 



■2^?'-'"-^'«')" = -?4^ "'''' 



Hq)lb{-q)\o) = ^E(oi«.(q)''-(-q)lo)-E.(q)-' 



•5^i:W".^.(q)-'^-Sj^ (B26, 
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(0|c(-q)la(q)|0) = ^ ^ ^(0|6,(-q)a.(q)|0)K(q)-' 



-^ Y. (rl„,^m^ciE^{q) ^ = - II 



a 






(B27) 



|o(q)^c(-q)|0) = ^^a(0|a.(q)6.(-q)|0)i?.(q)-i 



1 

'2S 



^o-/^,qm^qE^(q; 



-1 



BAq) 



^USE^iq)^ 



Then 



J 



12 



where 



W„ 



EE[^^(q)'e(q)]-^ 



(B28) 



'uc q o- 

[A - E.(q)][-7(q)e*^-'^ + 1 + 6(q)] + [A + E.(q)][-7(q)e^^-'^ + 1 - e(q)] 



-B^{q)[-j{q}e 

-5.(q)[-7(q)e 



-iqya/2 



-iqya/2 



+ (l + e(q))e 

:i - <q))e 



iqya/2-i iqa:a/2 



igya/2] iqxa/2 



-a5,(q)[-7(q)e^«^'^/2 + (1 + g(q))e-^'?«"/2]e*'?-'^/2 
-fffi,(q)[-7(q)e*«'''^/2 + (1 - e(q))e-^«^'^/2]e*«-"/2 
+a[A - i?.(q)][-7(q) + (1 + e(q))e^«-1 
+a[A + E4q)][-7(q) + (1 - e(q))e^^-'^] 



(B29) 



7iqj 



|[cos(g^a) + cos(gya)] 



(B30) 



and 



e(qr = l-7(qr. 



(B31) 



We use the fact that J3 ^ J. Only if a sum is divergent will it make a difference if we 
retain nonzero J3. So we tentatively assume no divergences and write 



B+{q) = 4:J cos{qxa/2) cos{qya/2) 
-B_(q) = — 4Jsin(g^a/2) sin(gj^a/2) . 



(B32a) 
(B32b) 
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We now simplify Eq. ( P29|) . We note that under the sum over wavevectors we can replace 
expi^iq^a) by 7(q). Let us apply the same reasoning to exp[2(g^ ± qy)a/2]: 



exp[i(ga. ± qy)a/2] = cos(ga.a/2) cos{qya/2) =F sin(ga,a/2) sin(gya/2) 

+i sm{qxa/2) cos{qya/2)a ± cos(ga;a/2) sin)gya/2) ] . 

After summation over wavevectors the imaginary parts will drop out. So 

exp[i{qx + qy)a/2] = cos{qxa/2) cos{qya/2) — sm{qxa/2) sm{qya/2) 

exp[i(g^ — qy)a/2] = cos(g^a/2) cos{qya/2) + sin(g^a/2) sm{qya/2) 



(B33) 



In this connection note that sums which are proportional to i?+(q)i?_(q) vanish. So 



(B34a) 



(B34b) 



^"■EE^ 1 



[A - E^iq)] [-7(q)' + 1 + e(q)] +[A + E^)] 



-7(q)' + 1 - e(q) 



+a [A - EM)] [7(q)e(q)] +cr[A + E^)] [-7(q)e(q)] 
- (^) [-7(q)^5.(q) + (1 + e(q))5.(q)] 



[-7(q)afi,(q) + (l-e(q))5.(q)] 



^^^) [-7(q)5.(q) + (1 + e(q))ai?.(q)] 
/^5<x(q)' 



4J 



-7(q)i?^(q) + (1 - e(q))ai?.(q)] 



^"EE^ 1 



8iv„,^^V^,(q)'e(q); V 

-J-i5.(q)^[l-a7(q)]) 



f2Ae(q)2 - 2K(q)e(q) - 2i^.(q)e(q)a7(q) 



(B35) 



Now we must understand how the wavevector sums are to be done. The unit cell is 



Sii = ax + ay , a2 = —ax + ay 



(B36) 



Thus the reciprocal lattice vectors are 
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Gi = {n/a){x + y) , G2 = in/a){-x + y) . 



Thus the sums are carried over the first zone, shown below in Fig. |^. 



2. Ao 



(B37) 



Thus 



A, 







Vl-ll — amGn 



0-0 



CtmCn — V/_// 







(B38) 



where we invoke perturbation theory relative to decoupled Cul and Cull subsystems. As for 
Ai effectively we have Eq. ( [B4| ). Thus 

1 



Then 



A2 = -Jl2SY.U) 

-JuSY^/o 



i^c 






-JuSYi^ 

ieb ^ 

-Ji25E(0 









^m r> ^i 



(^m p "i 



A, 



i^d 



J12S 



a -S 



0)(0|[eJ+,^, + /.+^.,]e,n+^„j0) 

0) mels,^ + fi+s,f]e„.+s JO) 
0) {0\e^+sjels.^ + f,+sj\0) 
0){0\e^+sjels,^ + f,+sj\0) 
0y0\em+sJels.^ + f^+5J\0) 
0) {0\em+sjel_^s,, + fi+s^^. 



a'[(l)-a[q} 



) e^^^^'^-''"*' 



x(0|[et(k)e*'^''- + /(-k)e*''''''/]e(k)e-'''-^-|0) 



J\2S 



EEio 



uc q,k ieb 



I 



&(-q)^a(q 



x(0|[enk)e 



t('^^p«k■<5b 



) e^('=i+^)''^i"^ 
/(-k)e*''''''-f]e(k)e-*^'''"=|0) 



(B39) 



52 






EE(o 



uc q,k iec 



c(-q)^a(q, 



q\ gi(q+k)-ri: 



X 



(0|[e^(k)e^^-''- + /(-k)e*''-''=-f]e(k)e-^''-^- 



■'^uc q,k i£d ^ ^ 



q\ g*(q+k)-ri 



x(0|[e1'(k)e*'^-'''*- + /(-k)e^^-'''^/]e(k)e-^^-'''"=|0) 



Jl2^ 



/V2 

uc q,k iGa 



EE(o 



«(q)^a (qj 



q\ g*(q+k)-ri 



x(0|e(k)e-^''-^"=[ef(k)e*''-''"^ + f{-k)e^^-^^f]\0) 



uc q,k iefc 



a(q)^&(-q) 



o\ e**^''^''^''"* 



x(0|e(k)e-*'^-^-[et(k)e*^-''''^ + /(-k)e*'^-^''^]|0) 
. 1 , 



uc q,k iec 



aiqj^c(-qj 



q\ g«(q+k)-ri 



x(0|e(k)e-*'^-^-[et(k)e*^-''- + /(-k)e*^-''^^]|0) 
^''^■EE(0a(q)>(q)0\e^('^+^)- 



uc q,k igd 



'^ 



x(0|e(k)e-^''-^"=[e"f(k)e^^-^'*= + /(-k)e*''-'''*-']|0) . 



Doing the sum over i we get 



A, 



J12S 



N, 



E(o 



uc 



J12S 

*uc q 



J12S 



E(o 

q ^ 



a (q)^«(qj 

^(-q)^a(q) 
c(-q)^a(q) 



^E(0rfnq)^a(q: 



*uc q 
J12S 



1 



■E(o 

*uc q 



iVu 






AT,, 



■E(o 

q ^ 

E(o 



«(q)^a'^(qy 
a(q)^&(-q' 



"■uc 



«(q)^c(-q) 



■f^E(Oa(q)>(q: 



'£ 



0\ (0|[e^(-q)e-*'i-''- + /(q)e-*'i-''''^]e(-q)e''i-''-|0) 
o\ (0|[e1"(-q)e~''i-'''"= + /(q)e-*''-''''^]e(-q)e*'i-''-|0) 
o\ (0|[et(-q)e-*'i-^- + /(q)e~^'i-^=^]e(-q)e''i-^ 

o\ (0|[e"f(-q)e-^'i-^''^ + /(q)e-^'i-''*]e(-q)e^'i-^-|0) 

o\ (0|e(-q)e*'i-^-[et(-q)e-*'i''- + /(q)e-*'i''" 
o\ (0|e(-q)e^'i-^-[e'f(-q)e-^'i''*"= + /(q)^^''-''' 
o\ (0|e(-q)e''i-''-[e^(-q)e-''i-^- + /(q)e-''*-^=/]|0) 

o\ (0|e(-q)e''i-''-[et(-q)e-''i-''''^ + /(q)e-''i-''*]|0) . 



This is 



(B40) 



(B41) 
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A,. 



J12S 



E(o 



«'iqj^«iqj 



'iVuc 

J12S 



E 



E 



•uc q 
JtoS 



^(-q)^a(q) 



1 



5:(0rft(q)-a(q) 

-''uc q \ ^ 



J12S 

J12S 



E 



1 



,t/ 



A^, 



E 



uc q 



J12S 



E(o 



'uc q 
Jl2»S 



aiqj^a iqj 



a(q)^&(-q) 



aiqj^c(-qj 



AT,,, 



5:(0a(q)-rft(q) 







]|et(-q)e(-q)|0) + (0|/(q)e(-q)|0)e'''-''/« 
(0|[et(-q)e(-q)|0)e*'i-''-'' + (0|/(q)e(-q)|0)e''i-(^--^^^) 
(0|[e^(-q)e(-q)|0)e''i-^- + (0|/(q)e(-q)|0)e*'i-^--*''-''^^ 
(0|[et(-q)e(-q)|0)e*'i-''-^ + (0|/(q)e(-q)|0)e*'i-(''--^^/) 
(0|e(-q)et(-q)|0) + (0|e(-q)/(q)|0)e^'^-^/«= 
"(0|e(-q)e^(-q)|0)e*'i-''-^ + (0|e(-q)/(q)|0)e^'i-(^--''''/) 
(0|e(-q)e^(-q)|0)e^'i-'''"- + (0|e(-q)/(q)|0)e*'i-(''--^-^) 
(0|e(-q)e^(-q)|0)e''i-^-''(0|e(-q)/(q)|0)e''i-(^-~''''/) 



(B42) 



Here the symbol 5fae denotes the vector which goes from an f site to an e site via an a site, 
such that fae is a sequence of nearest neighboring sites. So 

J12 



Ao 



4Ar„, 



EE f^^(q) 

q o- \ 



[A - EM)] [K - ^q^qe^^'^i + [^ + ^-(q)] 



fq fqvTtqC 



BM) 



^-iiqx+qy)a/2p _ ^ ^ ^i{qy-q^)a/2 



-a5,(q) y'^^y-^'^^'^/^ml - e-^('?^+'?'')'^/2/qmql - aB^) U^^y-^^y/^ll - (,-i{i^+iy)'^l^i 



.qmq 



-^ [A - E„] 



e *^""mq - /qmq 



a [A + ^,] 



e *^""mq - /qmq 



(B43) 



Making the same replacements as in Ai we get 
J 12 



A, 



AN 

^-' * uc q a 

-bm)bM){^JV 



J2 E ^-(q)"' ( [^ - Ea{q)] ml - /qmq7(q)J + [A + E^)] ll - /qmq7(q) 

5,(q)5,(q)(4J)-i 
-afi,(q)fi,(q)(4J)-i 



mq - 0-/qmq 



/q - 0-/qmq 



-afi,(q)fi^(q)(4J)-i 



am^ — IqiTiq^ 



al^ - l^m^ 



+a [A - E,] h{ci)m\ - l^m^ +a[A + E,] ^(q)/^ - l^m. 



(B44) 



This is 
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A, 



Ji 



12 



8iV„ 



q o- \ 



\A - -B,(q)l [l - f(q) - 7(q)'] +\A + EM)] [l + f(q) " l(l 



^'">'^\2-2.0<q)l-f;^^l|2.-2.(q)] 



V 4J 

+a[A-E^{ci)] 



So 



4J 



7(q)[l-e(q)]-7(q) 



+ a[A + E,(q)] 



7(q)[l + e(q)]+7(q) 



(B45) 



Ao 



Jl2 



EE 



8iv,,^^V^-(q)Mq) 



S.(q)^J-Ml-^7(q)] . 



2Ae(q)2 + 2E.(q)e(q) + 2E.(q)a7(q)e(q) 



(B46) 



3. Summary 



So 






where 



^""i^??G 



c - ^ 



EJqYe 



If we extend the sum over 

C -^ 



tie sum over —n/a < qx,qn < vr/a, then we 



4iV,, ^ \j 



2Ae{qf - B '-^^ ^- 



may write these as 

7(q)]) 



nfj-'[i- 



(B47a) 
(B47b) 

(B48a) 
(B48b) 

(B49a) 
(B49b) 



Of course, note that now X]q = '^Nuc- So it is convenient to introduce the notation ( )q to 
denote (2Ar„c)"^Eq- Then 
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Or 






p _ 1 / l-7(q)'-2co5'(a^^/2)coB>g,/2)|l-7(q)] \ 
4\ |1 - cos2(a5,/2) cos2(o,,),/2)] ^1 - 7(q)" /^ 

C,4( , ^+"<q> \ . ' (B51b) 



In the approximation that 7(q) = 0, etc. Ca = Cp = \. Numerical evaluation yields 

C„ = 0.1686, C^ = 0.4210. (B52) 

APPENDIX C: IN-PLANE GUI GUI INTERAGTION 

Here we reproduce by perturbation theory the gap found phenomenologically by Yildirim 
et al.y We treat an antiferromagnet on a square lattice (of lattice constant a) , in which there 
are two sublattices, a and h. The lattice is shown in the Fig. |^ with the magnetic unit cell 
within dashed lines. The magnetic unit cell has basis vectors 

ai = a^ + ar) , 

^2 = -ai + afj . (CI) 

We transform to bosons using Eq. (plSf) . 

First we consider terms Ti in the Hamiltonian which are quadratic in boson operators. We 
write 

n = nj + ns . (C2) 

Here 

Hj = 4JSj2(a\q)aiq) + 6t(q)6(q) + ^iq)[a\cip{-ci) + a(q)6(-q)]) , (C3) 
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with 



7(q) = ^[cos q^a + cos qyo] . 



(C4) 



and the sum over wave vectors is over the Brillouin zone associated with the magnetic unit celL 
Also 



ns = 6j,sJ2Ui^)-cy{k)] 



a(k)+a^(-k) 



6^(k) + 6(-k) 



(C5) 



where c^ik) = cos k^a and Cy(k) = coskya. 

Since the effect we wish to treat involves energies of relative order (l/S*), we now consider 
the fourth-order terms, V4, in the boson Hamiltonian, which we write as 



where 



V, = Vj + V,, 



{ij) 



'J J (^i ) 



(C6) 



(C7) 



where {ij) indicates that i is summed over a sites and j over nearest neighboring b sites and 



V5 = 6Ji Y^ as 

(ij) 



-iajajaj(6] + bj) — i(a| + ai)bjbjbj + olaMb^ 



ro 



(Q 



where a^ is +1 for x bonds and —1 for y bonds. 

We construct the effective quadratic Hamiltonian by taking all possible averages of pairs of 
operators out of the fourth order terms. Thus we have the effective quadratic terms 



(ij) 



aib]{{bj + 4)^) + 2ai{bj + a\){b]{bj + 4)) 
+2b]{bj + aj)((&j + a^i)ai) + {bj + a|)^(ai6]) 



(C9) 



and 



^Hs = SJi^as -|a|(6] + bj){a\ai) - \a\ai{a\{b] + bj)) - \a\a\{ai{b] + bj)) 

-\ai{b] + 6j)(a|aJ) - \{a\ + ai)bj{b]bj) - \b]bj{bj{a\ + Oj)) 
-\{a\ + ai)b]{bjbj) - \bjbj{{a\ + ai)b]) + a\ai{b\bj) + b]bj{a\ai) 
+a\bj{b]ai) + a\b]{aibj) + b'^jai{a\bj) + bjai{b^ja\) 



(CIO) 
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where (X) denotes an average with respect to the quadratic Hamiltonian. 

Since the quadratic Hamiltonian is real and Hermitian we can equate averages like (a|K) and 
{ttibj). Also at this order of (l/S) we only need keep Hermitian contributions to the effective 
Hamiltonian. Therefore we write 






(aib] + albj){{bj + alf) + 4(ai + b]){bj + al){b]{bj + a])) 



+{b] + aif{aib]) + {bj + 4)^(0^6]) 



(Cll) 



Next we consider Aif^. Here we can eliminate any terms which involve local averages (e. g. 
{a\ai)) because they multiply a function whose Fourier coefficient vanishes at zero wave vector. 
Thereby we have 



(ij) 



\alai{al{b] + bj)) - \a\a\{ai{b] + bj)) - \b]bj{bj{a\ + a,)) 
\bjbj{{a\ + ai)b]) + a\bj{b]ai) + a\b]{aibj) + b\ai{a\bj) + bjai{b]a\) 



(C12) 



Taking the Hermitian part of this we get 



(ij) 



-|a|ai(a|(6| + bj)) - |(aja| + aiai){ai{b] + bj)) - \b\bj(bj{a\ + Oj)) 



— \{b]b] + bjbj){{a\ + aj)6j) + {a\bj + aM){b\ai) + {a\b] + aibj){aibj) 



(C13) 



Thus we need the averages 

X2 



where 



{b]bj) = (aUi) 
Y,, = {b,al) = {b]a,) = Yo + a,Y 
Zij = {b^aj) = {bjai) = Zq + asZ , 

j j 



(C14a) 
(C14b) 
(C14c) 
(C14d) 



(C15a) 
(C15b) 
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J 3 



(C15c) 
(C15d) 



where the sums over j are restricted to sites that are nearest neighbors of site i. Now drop 
terms which sum to zero because of a^ and also those (hke X^ij <^5(Hbj) which do not contribute 
at zero wavevector. Then we get 



AHj + AH, 






(2Xi + AYo){aib] + ajbj) + Yo[ib]y + ^ + aj + [a] 
+4(X2 + Zo)(aJ + 6,)K + &l) 



i^^iE 



I] L 



(F + Z)[-Aa\a, - Ah]h, - a] - {a\f ~ b] - {h]f] 



+SY{a\hj + aih]) + SZ{a\h] + aihj] 



(C16) 



The coefficients can be evaluated straightforwardly. For instance, if one considers Tij as the 
unperturbed Hamiltonian and treats Tis as a perturbation, then one has 



j 

= E^5 



{()\h]a.i^n5\Q) + {ms\h]am 



(C17) 



where |0 > is the spin-wave vacuum and £ is the unperturbed energy of the virtual state. We 
give the evaluations 



Xi = AC2c{5.h/jf 
Fo = AC2d{U,/Jf 

Y = -8C2a{SJl/J) 

Z = -8C2bi6Ji/J) 



(C18a) 
(C18b) 
(C18c) 
(ClSd) 



where 



G 



2a 



1 

128iV 



E 



[ca;(q)-Cj,(q)]^ 



e(qj 



(C19a) 
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C2d = C2- C^c , (C19d) 

where e(q)2 = 1 - 7(q)2, C2 = Caa + C2b. 

To summarize: the effect of quantum fluctuations of the in-plane exchange anisotropy are 
contained in the effective Hamiltonian of Eq. (|C16| ). Since the result is given in real space, we 
can apply it now to the 2342 structure where it gives rise to contributions to the dynamical 
matrices written in Eq. (]69|). The terms proportional to X2 + Zq are taken into account by the 
spin-wave renormalization incorporated in Zc- 

APPENDIX D: IN PLANE ANISOTROPIC I II INTERACTION 

We start from Eq. ([7T|), which can be written as 

F12 = —6Ji2yS/2 2^ ejCj I aj+j, + a^^^ + di_x + di_^ — bi_y — bi^y ■ 

+5Ji2\J S/2 2_^fifi\ Ctj-x + O-i-x + di+x + C^j+x ~ ^i+y ~ h+y 



^i+y ^i+y 



^i—y ^i—y 



CjCjCi-i A s: 7 o / a /r, V^r rt i -C Ji Ji Ji-\ 



A6JuS^S/2j:[e] + e. - ^] - A5J^2S^ S/2Y.[fI + U- 



2S ' - V / ^L.* .» 2S 

i-yOi-y + C^j^yCiJ^y] 



+5Ji2\jS/2Y^{e\ + ei)[a\^^ai+x + 4-xC^J-x + ^l-yh-y + cj^„c,; 



+5Ji2\/S/2j2ifi + fi) [A-x(^t~x + 4+xC?i+x + b\+yh+y + c|-yCi-y] • (Dl) 

i 

Eliminate terms linear in the boson operators by the shifts 

°i *■ Cj -p S , Ji ^ Ji -\- S , 

ai ^ tti + t , hi ^ hi + t , Ci ^ Ci + t , di -^ di + t . (D2) 



The corresponding Fourier transforms are shifted by a factor y/N^. 



^ e, = e(0) -. e(0) + ^JN^cS , ^ a, = a(0) -. a(0) + v'iV,,t . (D3) 
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In what follows e will denote e(q = 0) and similarly for other operators. Then the linear terms 
in the Hamiltonian are 

V^ = -A6Ji2S^N^cS/2h + e + /t + /V (D4) 

The quadratic terms in the isotropic part of the Hamiltonian are 

V2 = (4J + 2J3)S{a^a + b% + c^c + Sd) 

+AJ2S{e^e + /V) + Ji2S\y + d^]f + [h^ + c^]e + [a + d]f + [h + cJeM 
+ (2J + 2.h)S{a^h'^ + c^S + ah + cd) + 2JS{a^c'^ + 6^rf^ + ac + hd) + 4J25(eV^ + e/) 
+ Jis^Ha^ + rf^e^ + P + c^]/^ + [a + d\e +[h + c]f\ . (D5) 

Now 

d{Vi + 1^2) 
(9e 



= \^NucS\ -A5Ji2^S/2 + 2Ji2t + 4J2S + 2Ji2t + 4J2S 



da 
'12 ^^ 4Jj2 



= jNucS (4 J + 2 J3)t + J12S + (2 J + 2J^)t + 2 Jt + J12S . (D6) 



For J% <^ 4 JJ2 we have 



.^^^^, (D7, 

0J2 



2Ji2g _ Jl2SJl2^JS/2 

8J + 4J3 J2(8J + 4J3) ■ ^ ^ 

As discussed in the text, these are the expected results. 

Now we record the terms in the Hamiltonian which are cubic in boson operators: 



7Y*^^) = 6Ju\/S/2l - Y^ ejci I aj+^ + ai+^ + d\_^ + di-^ - hi 



+2/ ^i—y ^i—y 



iee V / 

+ IZ(// + /i) ( (4-x(^i-x + 4+xdi+x + b^,+ybi+y + cVyCi-y j \ ■ (D9) 
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Now make the replacements of Eq. ( P3D to get 






~r Cti-^x ~r Oji—x ' ^i—x ^i—y ^i—y ^i+y ^i+y 



^i—x "•" '-^i-x + ^i-\-x ~r '^J+x i+y '^i+y ^i—y ^i—y 



+ Y.U} + fi){ 



ide 



'^i+y^i+y 



ige V 

iefV /J 



(DIO) 



Here we dropped the terms proportional to (a). They are smaller than those in (e) by 
</i2/(4J) ~ 1/50. Also, as before, to this order in 1/ S we may replace the perturbation by 
its Hermitian part. Then the effective quadratic terms above are 

+ O'i+x + "i-.T + di-x — 0,- 



H+x T ^i—x "T ^i-x ^i~y ^i—y '-^i+y Ci-^-y 



~r / \ Ji + Ji) I 0'^_^ + CLi—x + Wj_|_2. + CLi-\-x "i-i-»' "i-l-?/ C,_,, C,;_,, 
J6f V 



^j+j/ ^t+y ^i-y ^t-y 



+ U^l + (elr + 4eje,] + Y^iifD' + f^ + 4///.] 

'^ 2_^ I '^i+x'^i+a; + Cti-x^^i-x + Oj^_yOi^y + Cj_|_j^Q+y 1 
jge \ / 

+^ 2_^ I O-i-x^i-x + "j+a;"i+x + '^i+y'^i+y + '^i-y^i-y I f 






4J2 



E4 



(q)a(q) + 6t(q)6(q) + ct(q)c(q) + rft(q)rf(q) + et(q)e(q) + /^(q)/(q) 



+ 



+ 



(q)+6t(q)+ct(q)-dt(q) 



e(q)-/(q) + et(-q)-/t(-q) 



-a(q) + 6(q) + c(q) - rf(q) 



e^q) - /^q) + e(-q) - /(-q) 



+ et(q)et(-q) 



+ r(q)r(-q) , (Dii) 



which leads to Eq. (^. 

Now we look at the fourth order terms in the Cull - Cull isotropic exchange interaction. 
These are 



Vdm = -\J2 y\ elfJ.gfLsfi+s + fi+seUiCi + 2e\eif},J, 



2^2 ^ \ c-iJi^sJi+s 
i£e,S 



i+SJi+S 



(D12) 
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Substituting in two shifts of (e), this is 



VDM = -\{eYJ2Y. 



i,S 



4f.+s + 2flsUs + 24fls + iflsY + 2/. 



+5^i 



+fi+sel + 2eje, + e^ + 2eje, + 2flsh+5 + 2(eJ + e,)(/i., + /,+,) 



(D13) 



Taking the Hermitian part of this, we get 
{6Ji2?S 



V, 



DM 



I6J2 



E 

i,5 



e\h+s + ejls + 2e\fl, + 2e,/,+, + Afl^h+s 



+4ele. + \e1 + |(el)^ + \fl, + K/i,)^ + 2(el + e.)UU + /^+^) 



(5Jl2)'5 



I6J2 



E 



f 6e(q)/t(q) + 6et(q)/(q) + 8e(q)/(-q) + 8et(q)/t(-q) j (c. + c 

+ 16/Hq)/(q) + 16et(q)e(q) + 2et(q)et(-q) + 2e(q)e(-q) 
+2/t(q)/t(-q)+2/(q)/(-q 



(D14) 

which leads to Eq. (|75D. 

Contributions from quartic terms in the Cul - Cull interaction are smaller, i. e. of order, 
(5712)^^12/(^^2)5 if "^6 take out one factor of (e) and one factor of (a). Taking out two (a) 
factors gives an even smaller result. Taking out two (a) shifts from the Cul - Cul anharmonic 
term gives a contribution of order Ji2^Ji2/{.JJ2)- All these terms are neglected. 

APPENDIX E: IN - PLANE ANISOTROPIC II - II INTERACTION 

1. Self-Energy Due to Cubic Perturbations 

We start by discussing how one constructs the self-energy due to cubic perturbations. The 
point is that we wish to avoid the complexities involving Matsubara sums etc. Let us suppose 
that we have an unperturbed Hamiltonian in terms of normal mode operators, -E(q) and F{q)\ 



n = 5: ^(q) E\ci)E{ci) + Ft(q)F(q) 



(El) 



Now we want to identify the perturbative contributions to the matrices A(q) and B(q). Suppose 
we wish to calculate perturbative contributions leading to an effective quadratic Hamiltonian 
of the form 
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ii?(q)Et(q)Et(-q) . (E2) 

For this purpose we make the identification 

SBiq) = {0\E{q)Ei-ci)V^V\0) . (E3) 

Thus for a;(q) ^0 and considering only the ground state, we may write 

SB{ci) = {0\dV/dE\q)^dV/dE\~ci)\0) 

\{^\dVldE\-qy-dVldE\c^\^) . (E4) 

Similarly for the term in the Hamiltonian 

A(q)Et(q)F(q) (E5) 

we make the identification 

M(q) = (0|E(q)rivFt(q)|0). (E6) 

Thus for co'(q) -^ and considering only the ground state, we may write 

5A(q) = (0|9l^/5E^(q)^ai//aF(q)|0) 

+ (0|9y/9F(q)^9l^/9Et(q)|0) . (E7) 

This type of relation holds generally under the two assumptions: a) we consider the perturbation 
to modes whose energy can be neglected in the energy denominators, and b) we consider only 
the ground and low lying excited states, so that boson occupation numbers are zero. 

We have made the identification in terms of the normal mode operators, but equally we 
may transform to any set of modes. 

2. Application to Cull - Cull In-Plane Interactions 

We start from Eq. (^) and implement the results of the preceding subsection. For small k 
we write 
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T^ ^ dV/de{k) = 6jj2S/N^,Y.[fi^) + f\-^)]e\^){c. - Cy) 



T, ^ dV/de\\,) = 5J2v/25/iV,eE[/(q) + /n-q)]e(-q)C 



(^X Cy J 



T, ^ dV/df{k) = -6jj2S/N^,Y.Hq) + et(-q)]/t(q)(c. - Cy) 



n = dV/df\k) = -(5J2v/25/iV,eE[e(q) + en-q)]/(-q)(c. - Cy) , 



(E8) 



where Cx = cos{aqx) and Cy = cos(aqy). Thus if p denotes — p, then 

( Ti Iti) = ^-^%^ Y: ([/(q) + f^meHci){cx - Cy)hf{p) + /t(p)]et(p)( 



(^X Cy ^ 



2S{6J2)' ^ ^^^^^^ _ m^E^]i-m^F^)icx - Cy)^[-m^El + /pF^]/p4(c. - Cy] 



iV„, 



q,p 



2S{6J2f ^ {cx - Cy)%ml _ {6J,f ^/(c. - Cyf\l{qf 



AT,, 



E 



'uc q 8J2S'e(q) ^-'2JVuc 

where we introduced the normal mode operators 



El 



4J2iV„. VV e(q)3 



(E9) 



Ed = /qe(q) - m^p{q) , 

^q = ^q/(q) - "^qe"^(q) , 



(ElO) 



where 



2 1 + e(q) 



m^ 



1 - e(q) 



^ 2e(q) ' ^ 26(q) 

where 7(q) = |[cos(ag2.) + cos(agj,)]. Similarly 



/qmq 



7iqj 

26(q) 



(Ti^T^) = {T,^T^) 



4J2iVuc 



2; V^/ \Cx- Cy) 



E 



e(qj 



l-<q) I 7(q)^^ 



(Ell) 



(E12) 



1 



{dJ^f ^[{Cx-Cy)'\3^{qf 



(T'l-Ts) - (T3-T1) - -p— ^y— ^. 

£: £: 4J2A^uc q V e(q) 



(E13) 



^^4^^^^<4^^^ 47.iV, 



isj2r ^1 



(.Cx CyJ 



[1-^q)]^ 



'UC q V ^[(ir I 4 

Now we have the contribution to the coefficient of e^e, which we denote ^055, as 
^a55 = (Ti^T2) + (T2^Ti) = {bJljJ^ [-32C2a - 16C2fe] , 



(E14) 



(E15) 
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where C2a and C26 were defined in Eq. ( |C19|) . 

Likewise the contribution to the coefficient of e^ f which we denote 5a^Q is 

^ase = (ri^T4) + (r4^Ti) = [(^72)7^2] [32^2^ - I6C26] • 
Similarly, 5h^ is the contribution to the coefficient of ^e^e^, so that 



(E16) 



and likewise 



5655 = 2(Ti-Ti) = [{5J2f/J2\ [-16^2^,] 



5^56 = (Ti^Ts) + {tJ^T,) = [{5J2f/J2] [48C2,] 



(E17) 



(E18) 



APPENDIX F: INTERPLANAR ANISOTROPIC CUII - CUII INTERACTION 

1. Pseudodipolar Interactions 

In order to facilitate the evaluation of the lattice sums we parametrize the anisotropic 
exchange interactions bewteen the zth Cull spin in one plane and the nearest neighboring jth 
Cull spin in an adjacent layer. We introduce the indicator variable ai which is unity if i is 
on the e-sublattice and is —1 if i is on the /-sublattice. We also introduce a variable /ij to 
distinguish between the two nearest neighboring sites with the same value of ctj. Then for the 
interaction between nearest neighboring Cull spins i and j in adjacent CuO layers we use Fig. 
4 to write the principal axes as 



{ij) 



n 



f^^ 



1(1 + aiaj)fi - \{l - ai(7j)i 



Hijjjj 



n 



(ij) 



|(1 + aiaj)^ COS ip + |(1 — aiGj)-?! cos ip HiHj + zsm.ip 
\{1 + 0"i0"j),^sin?/' + |(1 — aiaj)fjsm.ip fiifij — zcosip 



(Fla) 
(Fib) 
(Flc) 



We also write 



Sj = —ai{S — alai)^ + J S/2{a\ + a.j)fj + iJS/2{ai — a\)zai , 
where, in this appendix, a, is the boson operator for spin i. Then we have 



(F2) 
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s. ■ hl^^ 



-\{S - a\ai)c{oi + aj)fiifij 



+ |(1 - aiaj)^i^jC\JS/2{a\ + a^) + isai^S/2{ 
-1(5 - a|ai)s(o-i + aj)fiifij 



ai - a : 



+ 1(1 - aiaj)fiiiJ,jS^JS/2{al + a^) - icai^js/2{a 



4) , 



(F3a) 



(F3b) 



(F3c) 



where c = cosip and s = sin-?/'. Then we have 



il-ij — "J 2^ -f^mi^i ■ "^rraj [^j ' TT'm] — "J /_^ -timJ-n 



(F4) 



m=l 



m=l 



where, at quadratic order, 

Ti = i(l - crj(Tj)(aJaj + a]aj) + ^(1 + aiaj){al + ai){aj + a^) 

T2 = -i(l + (Jiaj)<?{a\ai + aja^) + 1 1 ^(1 - o"iO"j)(o.l + ai)cfiifij + iais{ai - a|) 

X |(1 — aiaj){aj + aj)c^i^j + iajs{aj — Oj) 



(F5) 



and T3 is obtained from T2 by replacing sinip by — cos'?/' and cosip by sin-?/'. Thereby we get 
the site-diagonal contribution to the Hamiltonian as 



6n = AAKSj2alai , 

i 

where AK was defined in Eq. (^). 

The remaining contributions to the Hamiltonian are found from Eq. ( |F4|) to be 



(F6) 






+\K-,S 



|c^(l — aiaj){a\ + aj)(a] + a^) — aiajS^{ai — ai){aj — a]) 



+icsfiifij{(7j — ai){af + ai){aj — a]) 



+ 



(FT) 



where . . . indicates further terms in K^ obtained from those of K2 by replacing cos ip by sin ip 
and sin-?/' by —cosip and j is summed over Cull nearest neighbors in adjacent planes. For 
Qx = % = ^ the imaginary term gives zero contribution to the dynamical matrices. Then, the 
terms with a^ = aj give a contribution to the Hamiltonian of 
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i j:aj=cri L 



Ki{al + aj)(a] + a^) — K2S^{ai — a\){aj — a]) 



—K'^(?{ai — a\){aj — a]) 

The terms with a, = —aj give a contribution to the Hamiltonian of 

S 



6n 



4 



E E 



K2C^{ai + aj)(aj + a|) + K2S^{ai — a\){aj — a|) 



+7^35^(04 + ai)(aj + a]) + K3C^{ai — al){aj — Oj) 



(F8) 



(F9) 



The term in Eq. ( |i^'6| ) and the number conserving terms in Eq. (]^^ reproduce Eq. ( ^2a| ) and 
the other terms in Eq. ( p8|) reproduce Eq. (|82c| ). Equation (|F9|) reproduces Eqs. ( ^2b|) and 



(I82HD . 



2. Dipolar Interactions 

For the dipolar interactions it is convenient to construct the Hamiltonian explicitly rather 
than to identify it with the pseudodipolar interaction. We substitute Eq. ( p^ ) into the dipolar 
interaction to get 

T-f-ij = g^A^i^ [Si ■ Sj - 3(Si ■ iij){Sj ■ iij)] -^ -Sg^filr^j^iSi ■ rij){Sj ■ rij) 
2, ,2 



^9 f^B 



v 



-cri[S - alai]{^ ■ Vij) + JS/2{ai + al){fi ■ Vij) + ioiJS/2{ai - a\){z ■ r^j) 



-aj[S - a]aj]{i ■ fy) + yS/2{aj + a]){i) ■ Vij) + iaj^S/2{aj - a]){z ■ Vij) 



■ (FIO) 



Here we dropped the term in Sj • Sj which may be included in the isotropic Heisenberg Hamil- 
tonian. At quadratic order this gives 



^ = 2 E ^ii = E 



2r?. 



(Ticrj(a]aj + ajai)(^ • Vijf - |(ai + a|)(aj + a!j){fi ■ Vijf 



+\aiaj{ai - a\){aj - a]){z ■ Yijf - iaj{ai + a|)(aj - a]){z ■ Vij){fi ■ Yij) 



(Fll) 



We now consider what contributions this gives to the dynamical matrix for Qx = % = 0- Then 
the imaginary term can be dropped. For simplicity we truncate the sums to include only 
interactions between adjacent planes. Then we have 
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da55 - 2^ —3 — 

i6e tj 
'^«56 = 2^ —3 

^055 - 2^ —3 — 

jee ij 

(^Ose = 2^ —-3 



ii ■ f'ijf - l^^ziV ■ f'ijf - \Cz{z ■ fijf 



-\{ffh,f + \{z-hjf 



+ E^^^.(e-^^.r (F12a) 



iG/ ^3 



\{v ■ hjY + \{Z ■ hj) 



-\{v ■ hjf - \{z ■ f,jf 



(F12b) 
(F12c) 
(F12d) 



ie/ *^ 

where C2 = cos(g2c/2), i is a fixed site in the e sublattice, and the sum over j is restricted to 
the planes adjacent to site i. 

This interaction is neghgibly small except with respect to the lowest in-plane mode. So we 
only need the combination 



^(055 + &55 - ^56 - h&) = Y. 



3g'f,j,S 



jeir.z,j=±c/2 y 



(^jii-f'ijf -(^jiV-^ijfCz 



E 



3^'^i^ 



^jKp^ij ~r Vij) (^jCz\Xij Vij ) 



jeir.zij=c/2 «J 
Note that the sum over sites j in an adjacent plane from site i vanishes: 



(F13) 



E 



<JjX,ij 



E 



5 / J 5 

jeII:Zij=c/2 ^ij jeII:Zij=c/2 ^^j 



^^j _ Q 



(F14) 



Thus 



(^(055 + hb - am - he) = 6(1 + c^)g jj^bS 



o'jXij'i/ij 



jeir.Zij=c/2 ^ij 



(F15) 



APPENDIX G: INTERPLANAR ANISOTROPIC GUI - GUII INTERACTION 

For the Cul sites we introduce further indicator variables r (which tells the direction of the 
moment) and p (which discriminates between sublattices) such that r = p = 1 for an a site, 
— r = p = 1 for a 6 site, r = p = — 1 for a c site, and r = — p = 1 for a c? site. Then, from Fig. 
5, we have the principal axes for the sites i and j where i (j) is in the Cul (11) sublattice as 
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m 



Hi) ^ Pi'^j 
x/2 



7712 = ZCOSCJ) + 



i + nv 

piaj sin I 



m 



(ij) 



-z sin 6 + 



V2 
Piaj cos 



V2 






(Gla) 
(Gib) 
(Glc) 



In checking the above it is useful to note that changing the sign of either pi or aj induces a 
180° rotation about the 2;-axis. 

Also we use Eq. (|F2|) for the Cull spins and 



Si = Ti{S - a\ai)^ + \J S/2{ai + a|)i7 + iTi\J S/2{a\ - ai)z 
for the Cul spins. Thus if i labels a Cul spin we have 

-{S - a\ai) - \J S/2{ai + a\) 



mf^-) 


Si = 


PiTiCTj 

V2 


rhf^ 


S^ = 


1 

V2 




^f' 


s.= 


1 

V2 


- 



Pi(7js{S — altti) + piajsJ S/2{ai + a\) + iTiC\'S{a\ 



-piajc{S — a]aj) + piffjcJ S/2{ai + aj) — iTiSvS{a\ 
and if j labels a Cull spin we have 



— a,; 



.fe) 



m^ ■ Sj 



i) c _ 



{S - a^jaj) - Tiaj^JS/2{aj + a]) 



^fe) . S, = ^ 



PiTis{S — Ojaj) + piajsJ S/2{aj + a]) + i(7jCvS{aj 






PiTic{S — Ojaj) + piajcJ S/2{aj + a]) — iajSvS{aj — 



where c = cos (/) and s = sin (/). We now write 



3 3 3 

'^ij = Yl V^rn ■ Si] ^ [TJlra " Sj] = 5 ^ -^m^rr 
r?i=l m=l m=l 



and at quadratic order we have 

Ti = \Tiaj[a\ai + a]aj\ + \[ai + a\\[aj + a]] 
T2 = \Tiajs'^[a\ai + a]aj] 



+ 



^PiajS[ai + a-j -\ —{a!: - ai 



ITiC 

V2^ 



la^c , 



1 / t\ oO jC 

2pi(rjs{aj + a]) + -^ 



-I) 



(G2) 

(G3a) 
(G3b) 
(G3c) 

(G4a) 
(G4b) 
(G4c) 

(G5) 

(G6a) 
(G6b) 
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T'i = \Tia j(?[a\ai + a^a 



\piajc{ai + a|) j={a\ - ai) 



+ 



\piajc{aj + a] 



U ^^i^ 



/2 '^'^J ^3' 



We drop terms which do not contribute to the dynamical matrix for q^ 
find that 



(G6c) 



qy = Q and thereby 



ieiLjeii ieijeii I 



This result reproduces that of Eq. (^ 



(G7) 
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FIGURES 




c/2 



^ * 



^ X 



FIG. 1. Magnetic structure of 2342. The Cul spins (in sublattices abc, and d) are thick arrows 
and the Cull spins (in sublattices e and /) are thin arrows. The basis vectors for the magnetic unit 
cell are ai = a{x + y), a2 = a{x — y), and as = ^{ax + ay + cz). All spin directions are in the CuO 
{x-y) plane. The ^ axis is defined to be collinear with the spin directions. 
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FIG. 2. Nearest neighbor vectors connecting magnetic ions in a CuO plane. Cul spins are filled 
circles and Cull spins are open circles. Left: the vectors 5+ and (5_ between nearest neighboring Cul 
spins. Center: the vectors 6x and 6y which give the displacements of nearest neighboring Cul's relative 
to a Cull. Right: the vectors 25x and 26y between nearest neighboring Cull spins. 
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FIG. 3. Spin-wave spectrum for wave vector = q^c/(27r) along the c direction in the absence of 
anisotropy. Each mode is two-fold degenerate. The left-hand scale applies to the lower modes and the 
right-hand scale applies to the optical mode. Left: without spin-wave interactions. In this case one 
mode has zero energy for arbitrary wave vector in the c direction. Right: with spin-wave interactions. 
In the presence of easy plane anisotropy, the two-fold degeneracy is removed and only one mode 
(corresponding to rotation within the easy plane) is gapless at zero wave vector. When the four-fold 
in-plane anisotropy is also included there are no gapless modes. 
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^ f c 

FIG. 4. Interplanar CuII-CuII interactions. Left: a plaquette of Cull spins in one plane with a 
Cull neighbor in the adjacent plane over the center of the plaquette such that the isotropic CuII-CuII 
interaction is frustrated. Right: The principal axes for the exchange tensor of a spin in the e sublattice 



at O with a spin in the e sublattice at A. The directions of the axes are given in Eq. (8C). The axes 
for the interactions of the spin at A with other spins in the lower plane can be obtained by a rotation 
of coordinates. 
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FIG. 5. Interplanar CuI-CuII interactions. Left: a plaquette of Cul spins in one plane with a 
Cull neighbor in the adjacent plane below the center of the plaquette such that the isotropic CuI-CuII 
interaction is frustrated. Right: The principal axes for the exchange tensor of a spin in the e sublattice 
at O with a spin in the c sublattice at A. The directions of the axes are given in Eq. (|80|). The axes 
for the interactions of other pairs of CuI-CuII nearest neighbors in adjacent planes can be obtained 
by a rotation of coordinates. 
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FIG. 6. Unit cell of the square lattice. 
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FIG. 7. Brillouin Zone for the Square Lattice. 
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TABLES 



TABLE L Definitions of Parameters. Notation: 6 J = (J_l — J\\)/2. 



a 


Ca « 0.1686 


r 


c 


^ 


C2 ~ 0.01" 


t(") 
"^12 


(^ a J 12/ J 


Eq. (HD 


{SJi)VJ 


{6Jl2f/J2 


{5J2?/'h 


Eq. (Cl) 


Eqs. (H), (1) 


AK 


Ai^eff 


X = 7x 10-^1-3 










Eq. (H) 


Eq. (|7D 


Eq. fl) 











a) See ReL |TI 
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TABLE II. Estimated Values of Parameters from Experiment and Theory 



Parameter 


Values in meV 




From 


Experiment 


From Theory 




Value 


Reference"^ 


Value 


Reference" 


J 


130 ±5 





145 


m 


J3 


0.14 ±0.02 


m 






J12 


-10 ±2 


g TW 






■h 


10.5 ±0.5 


1 






AJi(T = OK) 


0.081 ±0.01 


TW 


0.04 


11,14 


AJi(T = 200K) 


0.068 ±0.011 









AJ12 






1.3'' 


m 


AJ2 


0.004 ± 0.004 


g, TW 


0.036 


E 


5Ji 


±0.04 


m, TW 


-0.02 


14 


SJ12 


±0.027 


1 


-0.015 


(b) 


6J2 






0.4 




Ai^eflF 






2.73 X 10"^ 


TW, Eq. (13) 


a 


0.135 


IS TW 


0.13 


TW, App. B 


T 


1.2 X 10"^ 


1|, TW 


~io-5 


(c) 


c 


7 X 10"^ 


1 


2.2 X 10"^ 


(c) 


e 






10-6 


(c) 



a) TW denotes this work. 

b) This is the contribution to 6J12 from dipolar interactions, which is much larger than that 
estimated from 6 J/ J ~ 1.5 x 10"^. 

c) Evaluation based on the relevant J's. 

d) Evaluated for the similar compound Ba2Cu304Cl2. 
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TABLE III. Renormalizations J -^ ZJ 



Quantity 
Renormalized to 

Refer to: 


J 

ZcJ 

(1 + 0.085/5) J 
Ref. |20 


^3 

Z^Js 

(1 - 0.2/5) J3 
Eq. {4C) 




AKesZi 

(1 - 0.2/5)i^eff 

Eq. (|99|) 


VJAJi 


ZgVJAJj 


(1- 0.2/5) VJAJj 
Ref. |24 



a) In the dynamics JJ3 -^ ZcZ^JJ^ = Z^JJ^, where we set Z^ = 0.77. 

b) In the dynamics J2Ai^eflf —^ ^c-^3<^2^-^eflf = Z^J2AKefi, where we set Z| = 0.77. 

TABLE IV. Experimental Values of Spin-wave Gaps at Zero Wave Vector 



Mode 


Temperature 


Energy (meV) 


Ref 


u.> 


T = 200K 


5.5(3)^ 


1 


L0> 


T = 200K 


0.066(4) 


16 


U> 


T^OK 


10.8(6) 


E 


UJ> 


T^OK 


9.1(3) 


Izl 


.< 


T^OK 


1.7473(4) 


16 


.< 


T^OK 


1.72(20) 


E 


UJ< 


T^OK 


0.149(3) 


16 



a) Extrapolated to T = 0. 
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TABLE V. Values in (10 ^meV) of the four-fold anisotropy constant k 



k 


T = 1.4K 


T = lOOK 


Experimental: From statical 

Experimental: Fitting Eq. (109|) to AFMR datall 

Theoretical: See Eq. (108|) 


25 

41 
56 


2 
1 

1" 



a) \J± — J\\\ = 0.041 meV is fixed so that the dynamics and theory agree. 

TABLE VI. Amplitude of the Dynamic Structure Factor 
Results for wave vector 2Ti{Hx/a + Ky/a + Lz/c) + g^ for H and K half integral and H + K + L 
an even integer. Results are given only to leading order in J. ys = x^, + AAJiS + |a = 
2J3S[l — cos{qzc/2)]+4:AJiS+^a. The mode energies (without 1/S corrections) and intensities 
[/"'^(q)] are independent of the particular values of H, K, and L and are evaluated for g^ = 0. 



Mode Energy 



Energy^ (meV) 



Intensity 



Formula 



Evaluation 



a;> = [8JSy3] 



1/2 



10.^ 



TZZ 
-*> + 













UJ, 



{8J2S 



4A J S + 2°(^^-^i-5'+^3) 



1/2 



1.72 



TZZ I6J2S 

III = 



12 




[8JS{x3 + 2a)] 



1/2 



9.1 



TZZ 













8J2S 



2ax3+Q4:a{2T-C)C2+?>aC,S 
xz+a 



1/2 



0.15 



TZZ 

jm _ 



16J25 
, .< 




140 
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a) See Table IV. 

TABLE VII. Amplitude of the Dynamic Structure Factor 
Results for wave vector 27T{Hx/a + Ky/a + Lz/c) + qz for H and K integers and H + K + L 
an even integer. The notation is as in Table V. Results are given only to leading order in J. 
The intensities are evaluated for g^ = and H = L = 1 and K = 0. 



Mode Energy*^ (meV) 



Formula for Intensity 



Intensity 



Lol = 10.8 



tVV 



TZZ 



8JS 



8JS 



1 



-1) 



y3[l + (-l)^] + (-l)^ 



a 



50 




col = 1.72 



TZZ 
-'< + 



tV^I 



' {8JSy{4J2S)a^ 
>4 < 

,< 



4J2S 



26 




u;> =9.1 



TZZ 



8JS 



l + (-l)^][x3 + ia] + «(-l)^}^ 
l-(-l)^ 



Tm 



8JS 




59 



uj< = 0.15 meV 



TZZ 






4J2 5 
(4J25)(8J5)^a^ 






1 





570 



a) See Table IV. 
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